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ABSTRACT 

The  present  work  is  an  investigation  into  the  scale  effects  on 
cavitation  Inception  in  submerged  water  jets.  Four  scale  effects  were 
studied:  (1)  jet  size,  (ii)  jet  velocity,  (iii)  dissolved  air  content, 
and  (iv)  the  nuclei  population  in  the  flow.  The  nuclei  population  was 
artificially  altered  by  "electrolysis  seeding." 

Holography  and  schlieren  photography  were  used  to  observe  the 
flow.  Direct  measurements  of  the  nuclei  population  were  also  accom¬ 
plished  by  holography.  In  addition,  the  instantaneous  pressure  field 
in  the  jet  was  successfully  mapped  out  using  specially  tailored 
bubbles  as  pressure  sensors. 

It  was  found  that  Inception  did  not  generally  occur  in  the  cores 
of  the  turbulent  eddies  and  that  the  region  in  which  the  cavities  were 
first  seen  were  dependent  on  the  size  of  the  jet. 

Pressure  measurements  showed  that  negative  peak  pressure 
fluctuation  intensities  of  as  high  as  120  percent  of  the  dynamic  head 
existed  in  the  jet.  The  results  also  revealed  that  the  Instantaneous 
pressure  fluctuations  have  a  slightly  skewed  bell  shape  probability 
distribution. 

For  the  present  tests,  the  inception  index  was  independent  of 
the  exit  velocity  for  a  constant  size  jet.  However,  when  the  flow  was 
seeded,  the  Inception  number  showed  a  linear  dependence  on  velocity 
and  this  dependence  Increased  with  Increased  number  of  seeded  nuclei. 
The  effects  of  the  nuclei  number  density  and  pressure  fluctuations  are 
Incorporated  in  a  "probable  cavitation  occurrence"  parameter  which 
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shows  promise  In  reducing  scatter  In  this  type  of  experiment. 


♦ 


» 


t 


# 


vi 


TABLE  OF  CONTENTS 

ACKNOWLEDGMENTS 

ABSTRACT 

LIST  OF  SYMBOLS 

LIST  OF  FIGURE  CAPTIONS 

CHAPTER  1:  INTRODUCTION 

1 . 1  General  Background 

1.2  Cavitation  Nuclei 

1.3  Scale  Effects 

1.4  Equipment  and  Instrumentation  that  are 
Currently  Employed  in  Cavitation  Research 

1.4.1  Flow  Visualization  and  Nuclei 
Measurements 

1.4.2  Detection  of  the  Inception  Point 

1.5  Review  of  Recent  Studies  on  Scale  Effects 

1.5.1  Flows  Around  Axi symmetric  Bodies 

1.5.2  Free  Shear  Flows 

1.6  Scope  of  the  Present  Work 

CHAPTER  2:  DESCRIPTION  OF  THE  EXPERIMENTAL  APPARATUS 
AND  INSTRUMENTATION 


Page 

ii 

iv 

ix 

xii 

1 

1 

2 

4 

6 

7 

10 

11 

11 

14 

16 

18 


2.1  General  Description  18 

2.2  Low  Turbulence  Water  Tunnel  [LTWT]  19 

2.3  Deaerator  System  19 

2.4  Flow  Apparatus  21 

2.5  Electrolysis  Bubble  Production  23 

2.6  Inception  Detection  System  24 

2.7  Holocamera  and  Reconstruction  System  24 

2.8  Schlieren  System 


27 


> 


ft 


ft 


ft 


t 


« 


♦ 


« 


vii 

TABLE  OF  CONTENTS  (continued) 

Figures  for  Chapter  2 

CHAPTER  3:  EXPERIMENTAL  PROCEDURES  AND  MEASUREMENT 
TECHNIQUES 

3.1  Introduction 

3.2  Preliminary  Preparations 

3.3  Measurement  of  the  Incipient  Cavitation 
Index 

3.4  Flow  Visualization  and  Nuclei  Measurements 
Figures  for  Chapter  3 

CHAPTER  4:  MEASUREMENT  OF  PRESSURE  FLUCTUATIONS  IN  JETS 

4.1  Introduction 

4.2  Principle  Behind  Static  Pressure 
Measurements  Using  Bubbles  as  Pressure 
Sensors 

4.3  Summary  of  Requirements 

4.4  Preliminary  Work 

4.5  Description  of  the  Bubble  Generator 
and  the  Injection  System 

4.6  Test  Procedures 

4.7  Data  Reduction 
Figures  for  Chapter  4 

CHAPTER  5:  PRESENTATION  OF  RESULTS 

5.1  Introduction 

5.2  Measured  Scale  Effects 

5.3  Flow  Visualization 


Page 

29 

36 

36 

36 

37 
39 

41 

42 
42 

44 

46 

47 

49 

50 

51 
57 
62 
62 
62 
66 


5.3.1  Observation  Under  Stroboscopic 
Lighting 


66 


Page 


r 

*  viii 

TABLE  OF  CONTENTS  (continued) 


5.3.2  Flow  Observation  with  Schlieren 
Photography 

69 

5.3.3  Flow  Visualization  with  Holography 

71 

► 

5.4  Nuclei  Population  Density  Measure¬ 
ments  by  Holography 

72 

5.5  Measurements  of  Pressure  Fluctuation 

75 

» 

Figures  for  Chapter  5 

78 

CHAPTER  6:  DISCUSSION  OF  RESULTS 

124 

6.1  Introduction 

124 

6.2  Unseeded  Jets 

124 

6.3  Seeded  Jets 

137 

6.4  Comparison  of  Experimental  Results 

141 

• 

6.4.1  Nuclei  Measurements 

141 

6.4.2  Photographic  Investigation  of  the 
Non-Cavitating  Jets -Coherent 
Structures 

142 

4 

6.4.3  Pressure  Fluctuation  Measurements 

144 

6.4.4  Cavitation  Results 

149 

6.5  Proposals  for  Future  Work 

151 

4 

Figures  for  Chapter  6 

153 

CHAPTER  7:  SUMMARY  AND  CONCLUSIONS 

162 

REFERENCES 

164 

t 

APPENDIX  1:  ORDER  OF  MAGNITUDE  ANALYSIS  OF  THE  TERMS  IN  THE 
RAYLEIGH-PLESSET  EQUATION 

174 

APPENDIX  2:  COMPUTATION  OF  THE  PROBABLE  CAVITATION 

OCCURRENCE  PARAMETER  ,  Np 

176 

• 

A2.1  Objective 

176 

A2.2  Illustrative  Example 

177 

i 

r 

4 


diameter  of  the  nozzle 
frequency  of  pressure  fluctuation 

volume  of  hydrogen  gas  produced  by  electrolysis  per  unit  time 

current  through  the  electrode 

constant  of  proportionality 

length  of  the  exposed  portion  of  the  electrode 

length  of  the  laminar  region  in  the  jet 

number  of  cavitation  events  per  unit  time 

total  number  of  nuclei  present  in  the  cavitation-prone  region 

of  the  jet 

number  of  nuclei  present  in  a  particular  size  category 
number  of  nuclei  that  are  normally  present  in  the  flow  field 
in  the  absence  of  flow  seeding 

number  of  hydrogen  bubbles  produced  by  electrolysis  per  unit 
time 

frequency  of  occurrence 
maximum  frequency  of  occurrence 
nuclei  number  density  distribution  function 
probable  cavitation  occurrence  parameter 
instantaneous  pressure  fluctuation 
cavity  pressure 

partial  pressure  of  the  gas  in  the  cavity  at  temperature  T 
reference  pressure  or  free  stream  pressure 
mean  static  pressure  in  the  water  tunnel  at  inception 
vapor  pressure  of  liquid  at  the  bulk  temperature 


X 


P„,P  (t)  instantaneous  local  static  pressure 

Xj  00 

P  partial  pressure  of  the  gas  in  the  bubble  at  the  reference 

9o 

position 

Pn  partial  pressure  of  the  gas  in  the  bubble  at  position  x 

cumulative  probabilities  of  favorable  negative  peak  pressure 
fluctuations  intensities  that  would  cause  the  explosive  growth 
of  a  nucleus 

probability  of  m  cavitation  events  occurring  per  unit  time 
rms  of  the  dimensionless  pressure  fluctuations 


P  (Pf) 


P  (m) 


Qj 

r 

R 

R_ 


m 


R 

R* 


R 

R 

S 


St 


volume  of  the  cavitation-prone  region 
volume  discharged  through  the  nozzle  per  unit  time 
radial  position  in  the  jet 
instantaneous  bubble  radius 

radius  of  the  bubble  immediately  after  exiting  from  the  nozzle 

radius  of  the  bubbles  produced  by  electrolysis 

geometric  mean  radius  of  the  ellipsoidal  bubble 

reference  bubble  radius 

radius  of  the  bubble  at  position  x 

mean  radius  of  the  bubble  or  nucleus 

D 

dimensionless  bubble  radius  (=  p-  ) 


jet  exit  Reynolds  number  (= 


V  .D 


) 


radial  velocity  of  the  bubble 

radial  acceleration  of  the  bubble 

surface  tension 

Strouhal  number  (*  t^-) 

j 


i 


t 

T 


oo 


or 

6 

Y 


v 

P 

CT 

Of 


0 


XT 


time 

bulk  temperature  of  the  liquid 
characteristic  bubble  response  time 
tangential  velocity  of  the  bubble  (=wQR0) 
exit  velocity  of  the  jet 
reference  velocity  or  free  stream  velocity 

c 

Weber  number  (=  *  „r-  ) 

°Rouo 

axial  position  in  the  jet 
dissolved  air  content 
average  dissolved  air  content 
damping  constant 
polytropic  index  (=  4/3) 

amplitude  of  oscillation  of  a  bubble  from  the  initial  size 

dimensionless  amplitude  of  oscillation  of  the  bubble  (=e/RQ) 

kinematic  viscosity  of  the  liquid 

density  of  the  water 

cavitation  number 

cavitation  inception  number 

circular  frequency  of  pressure  fluctuation 

natural  circular  frequency  of  the  bubble 

dimensionless  circular  frequency  (=  —  ) 


Fi gure 
Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


xii 

FIGURE  CAPTIONS 

2.1  Diagram  showing  the  layout  of  the  experimental  apparatus. 

2.2  (a)  Schematic  diagram  of  the  deaerator  system 
(b)  Photograph  of  the  deaeration  tank. 

2.3  Graph  displaying  the  performance  curve  of  the  deaeration 
system. 

2.4  Diagram  showing  the  geometrical  details  of  the  nozzle 
assembly. 

2.5  Diagram  illustrating  the  geometrical  shape  of  the  nozzle. 

2.6  Schematic  drawing  of  the  cavitation  inception  detection 
system. 

2.7  Diagram  of  the  holographic  reconstruction  system. 

2.8  Schematic  drawing  of  the  schlieren  system. 

3.1  Photograph  showing  the  mounted  nozzle  assembly  on  the 
top  window  of  the  water  tunnel. 

4.1  Diagram  of  the  air  bubble  injection  system, 

4.2  Diagram  showing  the  geometrical  details  of  the  bubble 
generator. 

4.3  Plot  of  e’0'^  cos  t  versus  t  (t^m 't). 

4.4  Photograph  of  the  reconstructed  image  of  a  string  of 


Injected  bubbles  in  the  potential  core  of  the  jet,  near 
the  nozzle  exit.  The  bubbles  are  128  urn  in  diameter. 


I 


» 


t 


♦ 


♦ 


xiii 

Figure  4.5  Photograph  of  the  reconstructed  image  of  ellipsoidal 
bubbles  in  the  transition  region  of  the  jet. 

Figure  4.6  Photograph  of  the  reconstructed  Image  of  a  badly  sheared 
bubble  in  the  transition  region  of  the  jet. 

Figure  5.1  Variation  of  the  incipient  cavitation  index  with  jet  exit 
Reynolds  number  for  the  3.17mm,  4,76  mm  and  6,35  mm  jet 
at  a  dissolved  air  content  of  14,1  ppm. 

Figure  5.2  Variation  of  the  incipient  cavitation  index  with  jet  exit 
Reynolds  number  for  the  3.17  mm,  4.76  irm  and  6.35  mm  jet 
at  a  dissolved  air  content  of  10,9  ppm. 

Figure  5.3  Variation  of  the  incipient  cavitation  index  with  jet  exit 
Reynolds  number  for  the  3.17  mm,  4.76  mm  and  6.35  mm  jet 
at  a  dissolved  air  content  of  7.7  ppm. 

Figure  5.4  Variation  of  the  incipient  cavitation  index  with  jet  exit 
Reynolds  number  for  the  3.17  nrn,  4.76  mm  and  6.35  mm  jet 
at  a  dissolved  air  content  of  4.3  ppm. 

Figure  5.5  Summary  of  the  variation  of  the  incipient  cavitation  index 
with  Reynolds  number  at  different  dissolved  air  contents 
for  the  3.17  mm  jet. 

Figure  5.6  Summary  of  the  variation  of  the  incipient  cavitation  index 
with  Reynolds  number  at  different  dissolved  air  contents 
for  the  4.76  mm  jet. 

Figure  5.7  Summary  of  the  variation  of  the  incipient  cavitation  index 
with  Reynolds  number  at  different  dissolved  air  contents 
for  the  6.35  mm  jet. 


xiv 


Figure  5.8  Variation  of  the  Incipient  cavitation  index  with  Reynolds 
number  at  different  electrolysis  currents  for  the  3.17  mm 
jet.  (a)  a  ■  13.98  ppm,  (b)  a  *  11.01  ppm,  (c)  a  ■  7.53 
ppm,  and  Id)  a  -  4.23  ppm. 

Figure  5.9  Variation  of  the  incipient  cavitation  index  with  Reynolds 
number  at  different  electrolysis  currents  for  the  4.76  mm 
jet.  (a)  a  ■  14.27  ppm,  (b)  a  =  10.79  ppm,  (c)  a  *  7.65 
ppm,  and  Id)  a  -  4.40  ppm. 

Figure  5.10  Variation  of  the  incipient  cavitation  index  with  Reynolds 
number  at  different  electrolysis  currents  for  the  6.35 
mm  jet.  (a)  a  =  14.0b  ppm,  (b)  a  =*  11.01  ppm,  lc)  7.92  ppm, 
and  (d)  3.97  ppm. 

Figure  5.11  A  plot  of  a.,  versus  ReQ  illustrating  the  effect  of  dissolved 
air  content  on  the  incipient  cavitation  index  when  the 
3.17  mm  jet  was  seeded,  (.a)  I  *  0.05  mA,  lb)  I  *  0.10  mA, 
lc)  I  *  0.20  mA,  and  id)  I  »  0.40  mA. 

Figure  5.12  A  plot  of  a.,  versus  Re^  illustrating  the  effect  of  dissolved 
air  content  on  the  incipient  cavitation  index  when  the 
3.17  mm  jet  was  seeded,  la)  I  =  0.6  mA,  (b)  I  3  1.0  mA, 
lc)  I  =  2.0  mA,  (d)  I  *  5.0  mA. 

Figure  5.13  A  plot  of  versus  Re^  illustrating  the  effect  of  dis¬ 
solved  air  content  on  the  incipient  cavitation  index  when 
the  4.76  mm  jet  was  seeded,  la)  I  =  0.05  mA,  (.b)  I  *  0.10 
mA,  (c)  I  a  0.20  mA,  and  Id)  I  ■  0.40  mA. 


XV 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


.14  A  plot  of  o^  versus  ReQ  illustrating  the  effect  of  dis¬ 
solved  air  content  on  the  incipient  cavitation  index 
when  the  4.76  nm  jet  was  seeded,  (.a)  I  =  0*6  mA, 

(bj  I  *  1.0  mA,  (_c)  T  3  2.0  mA,  and  (d)  I  -  5.0  mA. 

.15  A  plot  of  a.  versus  Re^  illustrating  the  effect  of  dis¬ 
solved  air  content  on  the  incipient  cavitation  index 
when  the  6.35  mm  jet  was  seeded,  (a)  I  *  0.05  mA,  (b) 

I  3  0.10  mA,  (c)  I  3  0.20  mA,  and  (d)  I  3  0.40  mA. 

.16  A  plot  of  versus  ReQ  illustrating  the  effect  of  dis¬ 
solved  air  content  on  the  incipient  cavitation  index 
when  the  6.35  mm  jet  was  seeded,  (a)  I  3  0.6  mA,  (b) 
1*1.0  mA,  (c)  I  3  2.0  mA,  and  (d)  i  *  5.0  mA. 

.17  Graph  of  a.,  versus  ReQ  showing  the  influence  of  jet 
size  on  the  incipient  cavitation  index  at  ct*14.1  ppm. 
la)  I  3  0.05  mA,  (b)  I  *  0.10  mA,  lc)  I  =  0.20  mA, 
and  Id)  I  *  0.40  mA. 

.18  Graph  of  versus  Re^  showing  the  influence  of  jet  size 
on  the  incipient  cavitation  index  at  a  *  14.1  ppm. 

(a)  I  *  0.6  mA,  (b)  I  *  1.0  mA,  (c;  I  *  2.0  mA,  and 
(d)  I  3  5.0  mA. 

.19  Graph  of  o^  versus  ReQ  showing  the  influence  of  jet 
size  on  the  Incipient  cavitation  index  at  a  *  10.9  ppm. 
(a)  I  3  0.05  mA,  (b)  I  3  0.10  mA,  (c)  I  *  0.20  mA, 
and  (d)  I  3  0.40  mA. 


Figure  5.20 

Figure  5.21 

Figure  5.22 

Figure  5.23 

Figure  5.24 

Figure  5.25 

Figure  5.26 
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Graph  of  versus  ReQ  showing  the  Influence  of  jet  size 
on  the  incipient  cavitation  Index  at  o  *  10.9  ppm. 

(a)  I  »  0.6  mA,  (b)  I  *  1.0  mA,  (c)  I  *  2.0  mA,  and 
Id)  I  *  5.0  mA. 

Graph  of  a.  versus  ReQ  showing  the  influence  of  jet  size 
on  the  incipient  cavitation  index  at  a  *  7.7  ppm. 
la)  1  a  0*05  mA,  lb)  I  a  0.10  mA,  (c)  I  *  0.20  mA, 
and  (d)  I  *  0.40  mA. 

Graph  of  a.  versus  ReQ  showing  the  influence  of  jet  size 
on  the  incipient  cavitation  index  at  a  *  7.7  ppm. 

(a)  I  *  0.6  mA,  lb)  I  *1.0  mA,  lc)  I  a  2.0  mA  and 
(d)  I  *  5.0  mA. 

Graph  of  versus  Reg  showing  the  influence  of  jet  size 
on  the  incipient  cavitation  index  at  a  *  4.2  ppm. 
la)  I  *  0.05  mA,  lb)  I  ■  0.10  mA,  (c)  I  »  0.20  mA, 
and  Id)  I  a  0.40  mA. 

Graph  of  versus  ReD  showing  the  influence  of  jet  size 
on  the  incipient  cavitation  index  at  a  =  4.2  ppm. 
la)  I  »  0.6  mA,  lb)  I  *  1.0  mA,  lc)  I  *  1.0  mA,  and 
(d)  I  =  5.0  mA. 

Schlieren  photograph  of  cavitation  inception  in  the 

C 

6.35  mm  diameter  jet  at  ReQ  ■  1.63x10  and  »  0.25. 
Schlieren  photograph  of  cavitation  inception  in  the  4.76 

5 

mm  diameter  jet  at  ReQ  ■  1.40x10  and  *  0.16. 
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Figure  5.27 

Figure  5.28 

Figure  5.29 

Figure  5.29 

Figure  5.30 

Figure  5.31 

Figure  5.32 

Figure  5.33 


Schlieren  photograph  of  cavitation  inception  in  the  3.17 

mn  diameter  jet  at  ReQ  *  1.02xl0S  and  *  0.15. 

Schlieren  photograph  taken  at  the  onset  of  the  almost 

'fcontinuous  popping  sound"  stage  of  cavity  development  in 

the  6.35  mm  jet  at  ReD  *  1.29xl05  and  a  *  0.19. 

U)  Schlieren  photograph  of  a  fully  cavitating  3.17  mm 

diameter  jet  at  a  ReQ  *  1.04x10  and  a  -  0.04. 

(bj  Schlieren  photograph  of  the  same  jet  as  in  (a)  above 

showing  the  region  where  the  cavity  collapsed. 

Schlieren  photographs  of  cavitation  inception  in  (a)  6.35 

mm  jet  at  ReQ  *  2.06xl0S  and  o.  *  0.22  lb)  4.76  mm  jet 
5 

at  ReQ  *  1.47x10  and  *■  0.09  lc)  3.17  mm  jet  at 

A 

ReQ  a  8.35x10  and  cn  *  0.08.  The  dissolved  air  content 
is  7.7  ppm. 

Schlieren  photographs  of  the  non-cavitating  6.35  mm  jet 
taken  with  the  knife  edge  parallel  to  the  flow  direction, 
(a)  ReQ  *  3, 03x1 04  Cb)  ReQ  =  7.20xl04  (c)  ReQ  = 

l.OOxlO5  Cd)  Re0  *  1.57x10s. 

Schlieren  photographs  of  the  non-cavitating  6.35  nro  jet 
taken  with  the  knife  edge  perpendicular  to  the  flow 
direction,  (a)  Re^  ■  1.06xl04  lb)  Re^  *  3.20xl04 
lc)  ReQ  »  6. 81x1 04  Id)  ReQ  =  1.64xl05. 

Schlieren  photographs  of  the  non-cavitating  4.76  mm  jet 
taken  with  the  knife  edge  parallel  to  the  flow  direction. 
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Figure  5 


Figure  5 


Figure  % 


Figure  5 

Figure  5 


Figure  5 

Figure  5 


(a)  ReQ  -  3.14X104  0>)  ReQ  =  5.94x104  (c)  ReQ 
»  8.2x1 04  (d)  ReQ  »  1.58x105. 

34  Schlleren  photographs  of  the  non-cavltatlng  4.76  mm  jet 
taken  with  the  knife  edge  perpendicular  to  the  flow  direc 
tlon.  (a)  ReQ  =  2. 1 3x1 04  lb)  ReQ  -  5.08xl04 

(c)  ReD  *  8.14xl04  Id)  ReQ  »  I.OlxlO5. 

35  Schlleren  photographs  of  the  non-cavitatlng  3.17  mn  jet 
taken  with  the  knife  edge  parallel  to  the  flow  direction. 

(a)  ReQ  »  1.50xl04  (b)  ReQ  =  3.57xl04 
(c)  ReQ  *  4.48xl04  (d)  ReQ  7.67xl04. 

36  Schlleren  photographs  of  the  non-cavitatlng  3.17  mm  jet 
taken  with  the  knife  edge  perpendicular  to  the  flow 
direction,  ta)  ReQ  *  1.41x1 04  (b)  ReQ  *  3.36xl04 

(c)  ReQ  *  7. 32x1 04  Cd)  ReQ  *  9.53xl04. 

37  Graphs  of  Lj/D  vs  ReQ  for  the  Ca)  6.35  mm  jet 

(b)  4.76  mm  jet  and  (c)  3.17  mm  jet. 

38  Plots  of  the  non-dimensional  jet  width  6/D  vs  the 
dimensionless  axial  position  for  the  (a)  3.17  irm  jet, 

(b)  4.76  nm  jet,  and  (c)  6.35  mm  jet. 

39  A  schematic  drawing  depicting  the  4  regions  in  the  jet 
where  nuclei  count  were  carried  out. 

40  Plots  of  the  nuclei  distributions  for  the  3.17  wn  jet. 

(a)  Region  1,  (b)  Region  2,  (c)  Region  3,  and  Id)  Region 
4. 
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Figure  5.41  Plots  of  the  nuclei  distributions  for  the  4.76  mm  jet. 

(a)  Region  1,  (.b)  Region  2,  (.c)  Region  3,  and  (d)  Region 
4. 


Figure  5.42  Plots  of  the  nuclei  distributions  for  the  6.35  mm  jet. 

(a)  Region  1,  (Jj)  Region  2,  (.c)  Region  3,  and  (.d)  Region 
4. 


Figure  5.43 


Figure  5.44 


Figure  5.45 


Figure  5.46 


Figure  5.47 


Figure  5,48 


Graphs  of  the  instantaneous  local  pressure  fluctuations 
as  a  function  of  the  dimensionless  axial  position  in  the 
3.17  mm  jet. 

Graphs  of  the  instantaneous  local  pressure  fluctuations 
as  a  function  of  the  dimensionless  axial  position  in  the 
4.76  mm  jet. 

Probability  histogram  showing  the  distribution  of  the 
pressure  fluctuation  intensities  in  a  jet  (Polytropic 
index  of  4/3). 

Probability  histogram  showing  the  distribution  of  the 
modulus  of  pressure  fluctuation  intensities  in  a  jet 
(Polytroplc  index  of  4/3). 

Probability  histogram  showing  the  distribution  of  the 
pressure  fluctuation  intensities  in  a  jet  (.Isothermal 
case). 

Probability  histogram  showing  the  distribution  of  the 
modulus  of  pressure  fluctuation  intensities  in  a  jet 
(.Isothermal  case). 
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Figure  6.1  A  plot  of  the  probable  cayitation  occurrence  parameter, 
Np.vs  a.  at  different  jet  exit  velocities  for  the  3.17 
mm  jet. 

Figure  6.2  A  plot  of  the  probable  cavitation  occurrence  parameter, 

Np,  vs  a.,  at  different  jet  exit  velocities  for  the  4.76 
ran  jet. 

Figure  6.3  A  plot  of  the  probable  cavitation  occurrence  parameter, 

Np,  vs  oi  at  different  jet  exit  velocities  for  the  6.35 
iron  jet. 

Figure  6.4  A  plot  of  Np  vs  for  the  3.17  mm,  4.76  ran  and  6.35  ran 
jet  at  different  jet  exit  velocities. 

Figure  6.5  A  plot  of  a.  vs  V.  for  the  equation  a.  *  K_  -  (8S/3pRV^). 

Figure  6.6  A  plot  of  the  probable  cavitation  occurrence  parameter, 

Np,  vs  cr..  at  a  jet  exit  velocity  of  20  m/s  for  the  three 
jets  when  flow  seeding  was  carried  out. 

Figure  6.7  A  plot  showing  the  nuclei  distribution  from  various 
sources  ^Katz  and  Acosta  (.1981)]. 

Figure  6.8  A  plot  showing  the  nuclei  distribution  in  the  present  jets 
at  different  dissolved  air  contents. 

Figure  6.9  Schlieren  photograph  showing  the  large  eddies  in  the  6.35 
mm  jet  at  ReQ  *  3.1x10  . 

Figure  6.10  Schlieren  photograph  showing  cavitation  inception  in  the 


core  of  a  vortex. 
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1.  INTRODUCTION 

1.1  General  Background 

One  of  the  first  questions  a  layman  or  a  newcomer  to  the  study 
of  cavitation  would  ask  would  be  "What  is  cavitation?"  Cavitation 
here  can  be  defined  as  the  occurrence  of  a  second  gaseous  phase  in 
a  liquid.  The  second  phase  could  be  a  mixture  of  some  noncondensable 
gas  and  the  liquid  vapor  or  just  pure  vapor.  The  change  of  phase  is 
usually  brought  about  by  dynamical  pressure  reduction  at  a  constant 
temperature. 

Among  the  first  persons  to  observe  cavitation  on  ship  propellers 
was  0.  Reynolds  in  the  late  nineteenth  century.  Since  then,  it  has 
captured  the  attention  of  a  wide  spectrum  of  professionals  ranging 
from  engineers  in  various  disciplines  to  applied  mathematicians, 
primarily  because  of  the  widespread  harmful  effects  that  accompany 
this  process.  The  performance  of  numerous  hydraulic  machines  has  been 
gravely  affected  by  cavitation,  especially  in  the  advanced  stages. 
Examples  are,  the  decreased  power  output  of  turbines  and  the  increase 
in  drag  on  marine  vehicles.  In  addition,  when  the  cavities  collapse, 
pressure  shock  waves  and  high  velocity  microjets  are  generated.  If 
these  shock  waves  and  microjets  strike  the  surfaces  of  machinery, 
erosion  by  pitting  of  these  surfaces  would  be  sustained.  Other  dele¬ 
terious  effects  of  cavitation  Include  noise  and  vibration.  As  a  direct 
result  of  all  these  adverse  effects,  cavitation  has  been  responsible 
for  very  expensive  premature  repairs  or  replacement  of  parts  of 
hydraulic  devices. 
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In  recent  years,  the  use  of  artificial  heart  valves  has  been 
on  the  increase.  The  valves  have  to  be  very  carefully  tested  to  Insure 
that  they  do  not  give  rise  to  cavitation  as  the  consequences  could  be 
fatal.  This  area  of  study  is  being  pursued,  for  example,  by  Freed 
etjl'39’. 

In  view  of  all  these  devastating  effects  that  accompany  the 
cavitation  process,  it  is  not  surprising  to  find  that  researchers  and 
engineers  alike  are  actively  engaged  in  extensive  related  research. 

It  is  hoped  that  with  increased  understanding  of  this  common  but  per¬ 
plexing  phenomenon,  better  future  designs  of  hydraulic  machines  and 
marine  vehicles  would  occur  with  eventual  huge  savings. 

A  first  step  towards  a  better  understanding  of  the  cavitation 
process  would  be  to  determine  the  cause.  This  would  require  a  close 
examination  of  the  microscopic  structure  of  the  liquid  itself  which  will 
now  be  briefly  considered. 

1 .2  Cavitation  Nuclei 

Plesset  has  shown,  from  his  theoretical  analysis,  that  a 
pure  homogeneous  liquid  could  sustain  tensile  stresses  of  several 
thousands  of  atmospheres.  However,  laboratory  tests  performed  by 
Dixon  (1909),  Vincent  (1943)  and  Scott  (1948)  as  reported  by  Knapp 
et  aj_  demonstrated  that  even  under  the  most  careful  tests, 
tensile  strengths  of  only  a  few  hundred  atmospheres  have  been  measured. 
This  discrepancy  in  the  measured  and  the  calculated  values  of  tensile 
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stresses  led  workers  in  this  field  to  postulate  the  existence  of 

"weak  spots"  in  real  life  liquids.  Two  questions  that  immediately 

come  to  mind  at  this  point  are  the  nature  of  these  "weak  spots" 

and  the  part  they  play  in  the  cavitation  inception  process. 

Solid  impurities  in  the  form  of  hydrophobic  particles  were 

proposed  by  Pease  and  Blinks  as  probable  "weak  spot"  candidates. 
(971  (541 

However,  Plesset  v  '  and  Hoi  1  v  1  have  examined  these  and  have 

dismissed  them  as  unlikely.  It  is  now  generally  accepted  that  these 

"weak  spots"  exist  in  the  form  of  undissolved  gas  or  vapor  pockets 

(341 

in  the  liquid.  Now,  Epstein  and  Plesset  v  '  have  demonstrated  that 
free  gas  bubbles  of  finite  sizes  would  rise  to  the  surface  of  the 
liquid  while  bubbles  that  were  small  enough  to  experience  Brownian 
motion  would  be  forced  into  solution  by  the  immense  surface  tension 
forces  acting  on  them.  Therefore,  for  the  undissolved  gas  or  un¬ 
condensed  vapor  pockets  to  exist,  some  agent  responsible  for  stabili¬ 
zing  the  bubbles  has  to  be  present.  Several  mechanisms  have  been 
proposed  to  explain  the  stabilized  states  of  these  pockets. 

Fox  and  Herzfeld  suggested  that  an  organic  skin  around 

the  gas  bubble  prevented  gas  from  diffusing  into  the  liquid.  Again, 

(971 

after  careful  consideration,  Plesset  u  1  has  concluded  that  this  was 

an  improbable  candidate.  The  most  widely  accepted  model  was  that 

(491 

suggested  by  Harvey  e£al_  u  ,  In  this  model,  Harvey  proposed  that 
the  undissolved  gas  could  be  trapped  in  the  crevices  of  suspended 
particles.  His  analysis  showed  that  when  the  liquid  pressure  drop¬ 
ped  below  the  vapor  pressure,  the  gas  pockets  trapped  in  the  crevices 
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would  "seed"  the  flow,  otherwise  they  remained  permanently  stabilized  in 
the  crevices.  A  further  analysis  of  Harvey's  model  and  a  critical 
review  of  other  proposed  models  can  be  found  in  Chapter  3  of  the 
monograph  by  Knapp  et  al_  and  also  the  work  of  Hoi  1 

With  the  nature  of  the  "weak  spots"  established,  the  remaining 
question  to  be  answered  is  the  role  they  play  in  bringing  about  the 
onset  of  cavitation  in  liquids.  It  is  obvious  that  since  the 
weakest  points  in  a  liquid  are  of  the  "weak  spots",  the  inception 
of  cavitation  would  take  place  at  one  or  more  of  these  spots.  Since 
these  spots  serve  as  the  nuclei  on  which  the  cavities  grow,  they  are 
commonly  referred  to  as  cavitation  nuclei.  Typically,  when  a  cavita¬ 
tion  nucleus  traverses  a  region  where  the  pressure  in  the  flow  is 
below  the  vapor  pressure,  the  bubble  will  explode  violently,  giving 
rise  to  the  term  "vaporous  cavitation".  It  is  conceivable  also 
to  think  of  gas  diffusing  into  the  bubble  across  the  interface  from 
a  super-saturated  liquid,  thus  giving  rise  to  a  "cavity".  This  type 
of  cavitation  is  often  referred  to  as  "gaseous  cavitation"  and  can 
take  place  at  liquid  pressures  both  above  and  be1^  the  pres¬ 

sure.  However,  this  process  proceeds  at  a  much  slower  rate  than 
"vaporous  cavitation".  Both  types  of  cavitation  have  been  observed 
experimentally. 

1 .3  Scale  Effects 

As  discussed  earlier,  cavitation  is  always  accompanied  by  many 
harmful  effects.  Consequently,  prior  to  the  construction  of  a  new 
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design  of  hydrodynamic  vehicle  or  hydraulic  machine,  it  is  desirable 
to  predict  its  performance.  If  the  design  is  known  to  cavitate  under 
normal  operating  conditions,  then  there  is  a  need  to  modify  the 
design  which  will  eventually  result  in  the  successful  operation  of 
the  machinery.  Unfortunately,  most  flows,  especially  those  that  are 
likely  to  cavitate,  are  very  complex  and  therefore  do  not  succumb 
readily  to  mathematical  analysis.  The  only  recourse  to  predict  the 
performance  of  the  new  design,  short  of  actually  constructing  a  pro¬ 
totype  and  operating  it,  is  to  do  model  tests.  These  tests  are  often 
performed  in  towing  tanks  or  water  tunnels  where  the  model  could  be 
mounted  and  tested  under  various  controlled  flow  conditions.  Once 
the  results  of  the  model  tests  are  obtained,  the  remaining  problem  is 
to  extrapolate  these  data  to  a  full  scale  model  or  prototype. 

An  important  parameter  in  cavitation  tests  is  the  cavitation 
number  defined  by 


where  and  are  some  reference  pressure  and  velocity,  respectively, 
P  ,  the  partial  pressure  of  the  vapor  phase  of  the  liquid  evaluated 
at  the  bulk  temperature  of  the  liquid  and  p,  the  density  of  the 
liquid.  If  the  tests  are  carried  on  in  a  water  tunnel  say,  the 
reference  pressure  and  velocity  are  usually  taken  to  be  the  free 
stream  pressure  and  velocity,  in  that  order.  The  cavitation  number 
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at  the  start  of  cavitation  is  referred  to  as  the  "incipient  cavitation 
index". 

It  follows  that  if  the  incipient  cavitation  index  is  the  only 
important  scaling  parameter,  then,  the  onset  of  cavitation  on  the 
test  model  and  the  prototype  would  occur  at  the  same  value  of  o  when 
both  are  operated  at  identical  dynamical  conditions.  However,  this 
is  found  not  to  be  the  case.  Cavitation  inception  has  been  found  to 
depend  not  only  on  the  model  size,  test  velocity  and  fluid  properties 
but  also  on  the  nuclei  population  in  the  fluid,  among  many  other 
parameters.  These  are  collectively  referred  to  as  "scale  effects". 
(Holl  and  Wislicenus 

Hence,  to  accurately  predict  the  performance  of  the  prototype 
from  model  tests,  a  thorough  knowledge  of  these  scaling  parameters 
and  the  corresponding  scaling  laws  is  required.  As  a  result,  a  lot 
of  attention  has  been  focused  on  this  particular  area  of  cavitation 
studies.  Two  types  of  flow  configurations  have  attracted  the  most 
attention  because  they  are  often  encountered  in  many  engineering 
situations.  These  are  (i)  the  flow  around  a  streamlined  body  of 
revolution  and  (ii)  free  shear  flows.  Recent  findings  made  in  these 
two  areas  will  be  discussed  in  Section  1.5  of  this  chapter.  For 
the  time  being,  the  equipment  and  instrumentation  that  are  presently 
available  for  cavitation  research  will  be  reviewed. 

1.4  Equipment  and  Instrumentation  that  Are  Currently  Employed  in 


Cavitation  Research 

Recent  advances  in  the  development  of  instrumentation  and  equip 


ment  used  in  cavitation  research  have  fallen  primarily  into  two 
categories.  The  first  category  is  the  methods  that  have  been 
developed  for  flow  visualization  and  nuclei  measurements.  The 
other  is  instrumentation  used  in  cavitation  detection.  These  will 
be  discussed  in  turn. 

1.4.1  Flow  Visualization  and  Nuclei  Measurements 

The  events  occurring  in  a  cavitation  process  normally  take  place 

very  rapidly,  and  this  makes  flow  visualization  difficult.  High  speed 

motion  pictures  have  been  successfully  employed  as  a  means  of  flow 

visualization  for  many  years.  An  example  is  the  now  classic  work  of 

Knapp  and  Hollander  in  recording  the  growth  and  collapse  of 

bubbles  on  an  axi symmetric  body.  Arakeri  and  Acosta  ^  succeeded 

in  using  the  schlieren  photography  to  look  at  laminar  flow  separation 

f45l 

on  the  hemispherical  nose  body.  Gates  v 1  later  used  the  same 

system  to  study  the  effects  of  transition  on  cavitation  inception  on 

several  axi symmetric  bodies.  The  "dual  photo  technique"  developed  by 
/  50} 

Hoyt  and  Taylor  v  1  was  utilized  to  look  at  cavitation  in  polymer 
jets.  Theirs  was  one  of  the  first  pictures  of  jet  cavitation. 

All  the  above  flow  visualization  techniques  provided  only  a 
two-dimensional  record  of  the  flow  field.  As  a  result,  much  important 
information  about  the  flow  was  lost.  This  problem  was  overcome  with 
the  advent  of  holography.  Among  those  who  have  employed  holography 
in  cavitation  research  were  van  der  Meulen  ^10)^  gates  et  al_  and 
Peterson  et  al_  Recently  Katz  use(j  a  combined  schlieren  and 
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holographic  technique  to  study  cavitation  associated  with  separated 
flows.  In  this  technique,  a  recording  of  a  slightly  heated  body  and 
flow  field  is  made  as  in  any  holographic  system.  The  difference 
enters  during  reconstruction  of  the  hologram.  A  knife  edge  is  used 
to  cut  off  the  main  beam  during  reconstruction,  thus  leaving  only 
the  diffracted  beam  to  be  observed.  In  this  way,  different  parts 
of  the  flow  could  be  analyzed  at  any  one  time,  by  correctly  positioning 
the  knife  edge,  without  sacrificing  the  rest  of  the  information.  This 
method  has  the  advantage  of  reproducing  a  good  image  because  of  little 
or  no  background  noise. 

It  is  known  that  nuclei  play  a  very  important  role  in  the 
inception  process.  However,  it  is  only  recently  that  they  are  becoming 
accessible  to  measurement.  The  size  and  density  of  these  nuclei 
have  been  measured  from  holograms  of  the  flow  field.  This  was 
carried  out  by  Peterson  et  al_  ^  and  Gates  ^45\  Unfortunately,  this 
is  a  very  time  consuming  process.  Another  method  developed  by  Keller^) 
to  measure  the  nuclei  population  involved  a  single  point  light 
scattering  technique.  The  principle  behind  the  operation  of  this 
system  is  to  measure  the  intensity  of  the  light  scattered  by  single 
nuclei  as  they  cross  a  focal  volume.  The  size  of  the  nucleus  is 
then  determined  from  an  earlier  calibration  curve.  This  method 
has  the  advantage  that  the  size  distribution  is  known  almost  in¬ 
stantaneously.  In  a  comparison  test,  Gates  et^  aT_  used  both  the 
holographic  and  a  light  scattering  system  which  is  similar  in 
many  respects  to  Keller's  ,  to  measure  the  nuclei 
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population  in  the  freestream  region  of  the  LTWT*.  They  found  that 
the  nuclei  population  distribution  measured  by  these  two  optical 
systems  do  not  agree  quantitatively.  One  reason  for  this  discrepancy 
is  the  inability  to  calibrate  the  light  scattering  system  accurately  with 
known  sizes  of  polystyrene  spheres.  Also  this  system  does  not  only 
register  a  count  of  the  bubble  size  it  sees,  but  also  records  a 
count  in  all  the  size  ranges  below  that  particular  size.  This  results 
in  a  nuclei  count  histogram  that  does  not  accurately  reflect  the  true 
size  distribution.  However,  studies  are  currently  being  carried  out 
at  ARL/PSU**  to  correct  the  short  comings  and  the  preliminary  results 
show  great  promise  as  reported  by  Yungkurth  and  Billet 
and  Kohler  and  Billet  ^72^.  Another  system  that  has  been  utilized 
to  measure  the  bubble  population  density  is  the  acoustic  doppler 
system.  This  system  is  described  in  some  detail  by  Oldenziel  ^86^. 

Comparative  tests,  such  as  the  much  reported  ITTC  tests, 

(Lindgren  and  Johnsson^8^)  have  revealed  a  poor  agreement  in  experi¬ 
mental  results.  This  apparent  inconsistency  is  not  only  in  the 
measured  inception  index  but  also  in  the  appearance  of  the  types  of 
cavitation.  Since  these  tests  were  carried  out  on  the  same  headform 
but  in  different  facilities  around  the  world,  it  was  concluded  that 
the  discrepancy  in  the  results  are  caused  by  "scale  effects".  In- 


*  The  Low  Turbulence  Water  Tunnel,  Keck  Laboratory  of  Hydraulic 
Research,  California  Institute  of  Technology, 

**  Applied  Research  Laboratory,  Pennsylvania  State  University. 
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disputed! y,  this  is  indeed  a  reason.  However,  this  lack  of  agreement  in 
the  experimental  results  could  also  be  partly  due  to  different  definitions 
of  the  inception  point  by  the  numerous  workers.  For  example,  one 
investigator  observing  the  flow  under  strobe-light  may  call  inception 
at  the  first  sign  of  a  cavity  while  another  may  interpret  it  different¬ 
ly.  Recently,  workers  in  this  area  are  becoming  more  aware  for  the 
need  of  a  universal  definition  for  the  inception  point.  To  accom¬ 
plish  this,  several  systems  have  been  developed  to  detect  the  onset 
of  cavitation  in  a  quantitative  manner. 

1.4.2  Detection  of  the  Inception  Point 

Two  systems  are  currently  employed  in  the  detection  of  the 
inception  point.  The  first  is  by  acoustical  means,  developed  by  a 
group  at  St.  Anthony  Falls  Hydraulic  Laboratory, and  the  other  is  by 
light  scattering.  Both  systems  involve  the  counting  of  cavitation 
events . 

In  the  acoustical  method,  a  pressure  transducer  is  mounted 
flush  on  a  test  body.  When  a  cavitation  event  occurs,  a  spike  is  ob¬ 
tained.  Inception  is  called  when  the  count  reaches  an  arbitrary 
predetermined  number  per  unit  time. 

Ellis  et  al_  and  Keller  used  a  light  scattering  technique 
to  detect  the  onset  of  cavitation.  In  this  case,  a  drop  in  the  light 
intensity  recorded  by  the  receiving  optics  was  taken  to  correspond 
to  a  cavitation  event.  Again  the  start  of  cavitation  was  denoted 
when  the  counted  events  reached  a  prechosen  value  per  unit  time. 
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Based  on  the  work  of  Ellis  and  Keller,  a  similar  system  was 
developed  and  tested  for  use  In  the  LTWT.  Just  as  In  Keller's  work, 
a  beam  of  light  from  a  He-Ne  laser  was  adjusted  to  graze  the  surface 
of  the  test  body  at  the  position  where  Inception  has  been  seen  to 
occur.  However,  the  beam  used  In  the  present  case  was  expanded  4 
times.  This  was  done  because  it  was  found  from  preliminary  work 
that  particles  and  bubbles,  not  associated  with  the  cavitation  incep¬ 
tion  process,  when  crossing  the  path  of  the  unexpanded  beam,  caused 
a  large  enough  drop  in  the  light  intensity  received  at  the  photo 
sensor  to  prematurely  trigger  a  count.  This  problem  was  overcome 
with  the  expanded  beam.  The  expanded  beam  was  eventually  focused 
onto  a  sensing  element. 

The  photo  sensor  used  was  a  "HP  PIN"  diode,  and  the  cavitation 
events  were  counted  on  a  frequency  counter.  A  count  of  50  events 
or  more  per  second  signaled  the  inception  point.  A  main  setback  of 
the  present  system  was  the  alignment  procedure.  Since  the  sensing 
area  of  the  PIN  diode  was  approximately  1  mm  in  diameter,  alignment 
was  critical  and  was  often  achieved  only  after  great  difficulty. 

Several  other  detection  systems  have  also  been  developed  by 
numerous  researchers  but  their  methods  are  not  as  widely  used  as 
those  described  above. 

1 .5  Review  of  Recent  Studies  of  Scale  Effects 
1.5.1  Flow  around  Axl symmetric  Bodies 

The  development  of  the  flow  visualization  and  Inception 


detection  equipment  described  above  has  served  as  an  added  impetus 
to  cavitation  research.  Recent  scale  studies  performed  with  the  aid 
of  these  instruments  will  be  briefly  reviewed.  It  must  be  emphasized 
that  the  following  discussion  is  not  exhaustive. 

In  the  quest  to  learn  more  about  the  scale  effects  on  cavita- 

(45) 

ting  flow  around  axi symmetric  bodies,  Arakeri  v  1  and  Gates  ^  1  used 
the  schlieren  system  to  study  the  effect  of  viscosity  on  flow  separa¬ 
tion  and  flow  transition  to  turbulence.  Their  work  demonstrated  the 
importance  of  viscous  effects  on  cavitation  inception. 

Peterson,  et  al_  and  Gates  and  Acosta  have  looked  into  the 
scaling  problem  associated  with  the  nuclei  population  of  the  flow. 

Their  holographic  studies  together  with  the  work  of  a  few  other 
investigators,  for  example,  Keller  who  used  a  light  scattering 
technique,  represent  the  relatively  few  papers  on  this  topic. 

Both  Gates  and  van  der  Meulen  ^ ^  looked  at  the  effect 
of  polymer  additives  on  the  cavitation  inception  process,  van  der 
Meulen  reported  finding  that  the  polymer  injected  appeared  to  sup¬ 
press  laminar  separation  on  the  hemispherical  body  by  initiating  a 
turbulent  transition.  This  resulted  in  a  delay  on  the  onset  of 
cavitation.  He  also  found  the  polymers  had  little  effect  on  the 
Schiebe  body;  a  body  which  does  not  exhibit  a  laminar  separation. 

These  findings  were  Independently  substantiated  by  Gates  ^). 

It  was  recently  suggested  that  conditions  on  a  prototype  could 
be  simulated  by  distributed  roughness  on  the  model.  This  prompted: 
several  investigations  on  the  effect  of  roughness  scaling.  Some 
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of  the  Investigators  who  have  undertaken  this  task  include  Billet  et  al 

Gates  et  al_  and  Oldenziel  The  group  of  investigators  at 

ARL  discovered  that  patch  cavitation  which  appeared  on  axisymmetric 
headforms  and  hydrofoils  could  be  controlled  by  surface  roughness. 

They  went  on  to  report  on  the  dependence  of  the  cavitation  index  on 
Reynolds  number.  In  general  it  was  found  that  the  cavitation  number 
for  a  given  body  size  increases  with  Reynolds  number,  and  increases 
with  decreasing  size  for  a  given  Reynolds  number.  Oldenziel  dis¬ 
covered  that  the  inception  number  on  a  Schiebe  body  with  attached 
roughness  was  unaffected  by  roughness  size  when  the  inception  point 
was  defined  as  the  instant  when  the  travelling  bubble  patch  reached 
a  length  of  one  cm.  This  is  indeed  quite  surprising.  Obviously, 
a  good  deal  of  work  has  still  to  be  performed  before  a  more  complete 
picture  could  be  acquired.  However,  the  results  reported  to  date  have 
provided  valuable  information  on  this  relatively  new  "scaling  effect". 

In  their  effort  to  explain  the  onset  of  cavitation  near  the 
surface  of  axisymmetric  bodies,  numerous  investigators  have  formulated 
theories,  for  example.  Parkin  Parkin  and  Kermeen  ^ ,  van  der 

Walle  Oshima  and  Johnson  and  Hsieh  These  and  other 

related  theories  will  not  be  discussed  here  as  they  are  beyond  the 
scope  of  the  present  work.  The  attention  will  now  be  focused  on  the 
studies  of  scale  effects  on  cavitation  inception  in  free  shear  flows, 
especially  in  submerged  jets,  which  is  the  prinary  objective  of  the 
present  Investigation. 

I 
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1.5.2  Free  Shear  Flows 

The  study  of  cavitation  in  free  shear  flows,  namely,  jets  and 
wakes,  is  of  great  interest  because  they  are  encountered  in  many 
engineering  applications,  such  as  in  the  separated  flow  around  pro¬ 
pellers  and  hydrofoils,  Shear  flows  are  also  found  just  downstream 
of  gate  valves  in  pipes  or  in  the  mixing  layers  that  are  formed 
downstream  of  control  gates  in  dams  (Amdt^^).  However,  unlike 
studies  of  flow  around  stream-lined  bodies,  relatively  little  is  known 
about  such  cavitating  flows. 

Kermeen  and  Parkin  carried  out  cavitation  tests  behind  sharp 
edged  disks  of  different  diameters.  They  reported  that  the  cavitation 
inception  number  was  a  unique  function  of  the  Reynolds  number.  However, 
Arndt  and  George  through  their  survey  of  numerous  literature, 

found  that  although  the  inception  number  increased  with  velocity  for  all 
the  model  sizes,  the  rate  of  increase  of  the  inception  index  with  velo¬ 
city  was  less  for  the  larger  disks  than  for  the  smaller  ones. 

The  first  person  to  investigate  cavitation  inception  in  jets 
was  Rouse  ^3),  sjnce  then,  numerous  other  investigators  have 
become  seriously  involved  in  the  studies.  The  following  table  is  a 
partial  listing  of  the  investigators  and  the  studies  undertaken  by  them. 

Name  of  Investigator 

Arndt  and 
Arndt  et  al 
C8, 9, 1571271 3) 


_ Scale  Effects  Studied _ 

Air  content,  jet  size,  polymer  additives 
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Baker  et  al 
041 

Ball  05} 

Hoyt  C58} 

Jorgensen  C631 

Lienhard  and 
Stephenson  (78} 


Atr  content,  polymer  additives 

Jet  size 
Polymer  additive 
Jet  size 

Jet  size,  water  temperature 


Table  1.1.  A  tabulation  of  the  investigators  and  the  studies  per¬ 
formed  by  tffi^OTTaVTCaTi  ng  jet?: - 

A  quick  review  of  the  findings  of  these  researchers  showed  that  the 
measured  value  of  the  cavitation  number  did  not  seem  to  exhibit  an 
upper  limit.  Furthermore,  the  incipient  index  for  the  same  jet  size, 
differs  considerably  from  investigator  to  investigator.  To  date,  no 
one  has  succeeded  in  finding  a  universal  scaling  parameter  although 
a  few,  for  example,  Lienhard  and  Stephenson  and  Arndt  ^12^,  have 
managed  to  obtain  semi -empirical  relationships  for  some  of  the  inception 
data. 

A  possible  reason  for  the  differing  results  could  be  the 
different  nuclei  population  in  the  test  flows.  Few  workers  have 
looked  into  this  aspect  of  the  flow.  In  addition,  many  workers  have 
failed  to  document  carefully  the  region  in  the  jet  where  cavities 
were  first  seen,  for  example,  whether  inception  was  in  the  shear 
layer  between  the  potential  core  and  the  jet  boundary  or  in  the  region 
beyond  the  potential  core  where  the  shear  layers  have  merged.  It  is 
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possible  that  the  discrepancy  in  the  measured  inception  index  could  be 
due  to  the  fact  that  the  onset  of  cavitation  occurred  in  different 
regimes  of  the  jet  which  are  governed  By  different  scaling  laws.  It 
is  also  equally  Important  to  look  into  certain  features  of  the  flow, 
such  as  the  exit  condition  of  the  jet  and  the  spread  rate,  as  these 
may  provide  insights  into  the  initial  vortex  roll-up  and  the  size 
of  the  coherent  structures,  which  could  in  turn  affect  the  cavitation 
inception  process. 

In  brief,  in  order  to  achieve  an  understanding  of  the  scale 
effects  and  the  scaling  laws  that  govern  cavitation  in  shear  flows, 
one  must  not  only  carefully  study  the  environment  but  must  also  look 
into  the  fluid  dynamics  of  the  flow. 

1 .6  Scope  of  the  Present  Work 

The  points  raised  above  are  only  some  of  the  questions  that  need 
to  be  considered  and  answered.  It  was  with  these  thoughts  in  mind, 
that  the  present  study  was  undertaken.  The  aims  of  the  present  work 
were  two-fold.  The  first  was  to  measure  the  cavitation  inception 
number  of  jet  flows  under  widely  different  flow  conditions.  Four 
parameters  were  considered.  They  were  Ci)  jet  size,  (.ii)  dissolved 
air  content,  (ili)  nuclei  population  and  Ci v 1  jet  velocity.  The 
second  aim,  of  course,  is  to  draw  these  observations  together  in 
some  meaningful  way. 

Flow  visualizations  were  made  with  schlieren  photography  as 
well  as  holography.  Holography  was  also  used  to  determine  the  nuclei 
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population  in  the  various  regions  of  the  jet,  A  method  was  devised 
to  measure  the  pressure  fluctuations  in  the  jet.  This  method  in  brief, 
uses  wel 1 -characterized  air  hubbies  as  pressure  sensors. 

The  results  of  all  these  measurements  will  be  presented  and 
discussed  in  later  chapters  together  with  an  interpretation. 
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2,  DESCRIPTION  OF  THE  EXPERIMENTAL  APPARATUS  AND  INSTRUMENTATION 


2,1  General  Description 

•  The  experimental  objectives  required  an  experimental  setup 
wherein  a  jet  of  water  could  Be  discharged  and  made  to  cavitate  under 
prescribed  conditions.  This  not  only  called  for  a  reservoir  of 

•  water  which  could  be  treated  to  meet  the  test  requirements,  for 

example,  deaerating  the  water  to  the  desired  air  content,  but  also 

★ 

an  environment  where  the  flow  could  be  observed.  The  LTViT  was 
chosen  as  the  experimental  setup  because  it  met  the  above  requirements 
In  the  present  work,  the  jet  was  created  by  discharging  tunnel 
water  itself  from  a  pressure  vessel  through  a  contoured  nozzle  into 
the  center  line  of  the  LTWT.  Cavitation  inception  in  the  jet  was 
brought  about  by  lowering  the  pressure  in  the  tunnel  while  main¬ 
taining  a  constant  jet  exit  velocity.  Onset  of  incipient  cavitation 
in  the  jet  was  detected  by  a  microphone  placed  against  the  top  window 
of  the  tunnel.  Two  optical  systems  were  used  for  flow  visualization. 
These  were  the  holocamera  and  the  schlieren  system. 

A  general  view  of  the  experimental  apparatus  is  shown  in 
Figures  2.1.  To  facilitate  presentation,  the  various 
components  that  make  up  the  apparatus  are  discussed  separately, 
commencing  with  the  familiar  LTWT. 


*  see  Section  2.2  to  follow 
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2.2  Low  Turbulence  Water  Tunnel  Q.TWT1 

Ttls  facility  was  described  In  great  detail  by  Gates  ^5^-,  and  it 
will  be  redundant  to  repeat  that  description  here,  For  the  present 
work,  it  is  sufficient  to  say  that  the  pressure  in  the  tunnel  could 
be  lowered  by  the  action  of  a  vacuizn  pump.  The  side  windows  of  the 
8  1/3  ft  test  section  are  made  of  ground  and  polished  plate  glass 
thus  facilitating  flow  visualization.  The  tunnel  is  also  equipped 
with  a  filter  system  and  this  permits  the  water  in  the  tunnel  to  be 
kept  clean. 

One  of  the  aims  of  the  present  work  was  to  study  the  effects  of 
air  content  in  the  water  on  cavitation  inception.  This  required 
running  tests  at  different  air  content,  from  nearly  saturated  water 
to  water  with  very  low  air  content.  When  Gates  performed  his  tests, 
the  method  he  adopted  in  reducing  the  air  content  was  to  cavitate  the 
tunnel,  stop  it,  then  bleed  the  lines.  This  was  repeated  over  a 
period  of  many  hours.  He  reported  achieving  an  air  content  of  7  ppm 
after  operating  for  16  hours  in  this  manner.  This  method  was  not 
only  undesirably  time  consuming,  but  it  was  difficult  to  attain  an 
air  content  of  much  less  than  7  ppm.  For  this  reason,  a  deaerator 
system  was  designed  and  installed  recently  as  part  of  the  LTWT.  The 
following  is  a  detail  description  of  this  unit, 

2.3  Deaerator  System 

Figure  2.2(a)  is  a  schematic  diagram  of  the  deaerator  system, 
while  Figure  2.2(b)  shows  a  photograph  of  the  deaeration  tank.  The 
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main  component  is  a  closed  cylindrical  yessel  measuring  2,54  m 
ClOO"l  in  length  and  0,3.1  m  (3  ft  1  in  diameter.  Water  is  pumped 
(by  the  supply  pump)  to  the  top  of  the  tank  from  the  water  tunnel.  The 
water  is  then  forced  through  a  set  of  10  "Bete  Spiral  Nozzles",  Model 
#TF28FC,  mounted  on  a  horizontal  section  of  the  pipe,  suspended 
inside  the  tank,  from  the  top.  These  commercially  designed  nozzles 
effectively  atomize  the  water.  The  fine  sprays  then  rain  down  on  a 
packed  bed  of  "Plastic  Super  Intalo.x"  tower  packings  resting  on  a 
steel  grating  stretched  across  the  entire  length  and  width  of  the  tank. 
This  is  aimed  at  increasing  the  surface  area  from  which  gas  diffusion 
can  take  place.  The  water  flows  through  the  packed  bed  and  collects 
at  the  base  of  the  tank  where  a  second  pump  returns  the  water  back  to 
the  tunnel.  The  deaeration  process  begins  when  the  pressure  inside  the 
vessel  is  reduced  by  a  vacuum  pump. 

During  deaeration,  the  tunnel  water  is  continuously  circulated 
through  the  deaeration  vessel.  It  is  desirable  to  prevent  flooding 
of  the  packed  beds  because  this  will  drastically  reduce  the  surface 
area  for  gas  diffusion.  Therefore,  there  is  a  need  to  be  able  to 
adjust  the  water  level  in  the  vessel.  This  is  accomplished  by 
throttling  the  flow  on  the  supply  and  drain  lines  to  the  tank  until 
the  desired  water  level  is  achieved. 

To  prevent  the  drain  pump  from  running  dry,  the  drain  pump 
automatically  kicks  off  when  the  water  level  in  the  tank  becomes  too 
low.  A  check  valve  in  the  drain  line  stops  any  back  flow  into  the 
tank.  This  pump  remains  shut  off  even  though  the  supply  pump  is 


still  operating,  until  tfie  water  leyel  rises  to  a  predetermined  level. 
By  the  same  token,  the  supply  pump  shuts  off  when  the  water  level  in 
the  tank  Becomes  too  high  or  overflooding  will  occur.  The  supply 
pump  comes  on  again  when  the  water  drops  to  a  safe  level.  These 
"ON-OFF"  actions  of  the  two  pumps  are  controlled  by  a  4-level 
"Gems  Level  Switch",  Model  ILS800,  mounted  on  a  side  viewing  glass, 
and  a  set  of  latched  relays.  When  the  water  level  in  the  tank  is 
properly  adjusted,  both  pumps  will  run  continuously  throughout  the 
deaeration  process. 

The  present  deaeration  system  can  circulate  200  gpm.  It  has 
been  found  that  at  this  flow  rate,  the  deaeration  time  was  reduced 
when  the  water  tunnel  was  kept  at  a  partial  vacuum.  Figure  2.3 
shows  the  performance  curve  of  the  present  system  at  two  different 
tunnel  pressures.  It  should  he  pointed  out  that  the  lower  curve 
represents  the  limit  of  the  present  system  because  further  reduction 
of  the  tunnel  pressure  would  soon  cause  the  supply  pump  to  cavitate. 
The  air  content  is  measured  with  a  Van  Slyke  Blood  Analyzer. 

The  present  deaerator  thus  allows  an  easy  control  of  the  air 
content  in  the  tunnel  water, 

2.4  Flow  Apparatus 

The  flow  apparatus  is  made  up  of  a  cylindrical  pressure  vessel 
and  a  nozzle  assembly.  This  is  shown  schematically  in  Figure  2.1  . 

The  pressure  vessel  is  0.3  m  in  height  and  approximately  1.22  m  in 
width.  A  small  pump  is  used  to  pump  water  from  the  water  tunnel  to 
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this  vessel,  and  the  tank  is  pressurized  with  the  Institute  compressed 
air  supply.  The  inlet  pressure  is  regulated  with  a  pressure  regulator 
and  can  be  varied  from  5  psig  to  85  psig.  If  desired,  the  water  in 
the  pressure  vessel  can  be  heated.  This  is  done  by  circulating  the 
water  through  a  heat  exchanger  (which  consists  of  a  coil  of  copper 
tubing  sitting  in  a  hot  water  bath.)  which  is  located  upstream  of  the 
vessel.  The  flow  rate  through  the  nozzle  is  measured  by  a  "Flotech" 
turbine  flow  meter  which  is  placed  on  the  discharge  side  of  the 
pressure  tank.  The  pressure  vessel  is  connected  to  the  nozzle  assembly 
which  is  mounted  on  the  top  window  on  the  tunnel  by  a  1.9  cm  (0.75  inch) 
I.D.  pressure  hose.  A  ball  valve  in  the  line  allows  the  flow  in  the 
jet  to  be  shut  off  in  between  runs. 

The  nozzle  assembly  is  made  up  of  five  different  sections. 

They  are  (i)  a  0.25  m  (9")  straight  vertical  section  of  pipe,  (ii) 
a  vaned  long  radius  elbow  to  assist  in  turning  the  flow  through  a 
right  angle,  (i i i )  a  honeycombed  section  of  pipe  downstream  of  the 
elbow,  (iv)  a  settling  chamber,  and  (v)  a  circular  nozzle.  The 
geometrical  details  of  the  nozzle  assembly  are  shown  in  Figure  2.4. 

These  various  parts  will  now  be  described. 

The  elbow  contains  2  vanes  made  up  of  cut  out  brass  sheets  of 
0.76  mm  (0.03")  thickness.  They  are  equally  spaced  and  are  held  in  place 
by  wedging  them  in  the  inner  walls  of  the  elbow.  The  elbow  is  welded 
at  one  end  to  the  short  vertical  section  of  pipe  and  the  other  end  to 
the  honeycombed  section.  The  honeycomb  sits  in  a  shouldered  region  and  is 
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prevented  from  being  washed  downstream  when  the  settling  chamber  is 
in  position.  The  settling  chamber  could  be  unscrewed  from  the  honey¬ 
combed  section,  thus  revealing  the  honeycomb.  This  facilitates  the 
cleaning  and  replacement  of  the  honeycomb.  Finally,  nozzles  of  dif¬ 
ferent  sizes  can  be  attached  to  complete  the  assembly.  Three  dif¬ 
ferent  nozzle  sizes  were  used.  They  were  3.17.  mm  (0.125  inch),  4.76 
mm  (0.188  inch)  and  6.35  mm  (0.25  inch)  in  diameter.  The  shape  of  a 
typical  nozzle  used  in  the  present  work  is  illustrated  in  Fig.  2.5. 

2.5  Electrolysis  Bubble  Production 

Electrolysis  bubbles  have  been  used  in  the  past  by  numerous 
investigators  for  flow  visualization.  In  the  present  studies, 
electrolysis  was  used  to  alter  artificially  the  nuclei  population  in 
the  jet  to  levels  which  are  well  above  what  is  normally  found  in  water. 

A  1.59  iron  O.D.  tungsten  rod  was  used  as  the  electrode.  It  was  com¬ 
pletely  insulated  with  "heat  shrink"  tubing  with  an  exception  of 
about  0.5  mn  of  the  pointed  tip.  The  electrode  was  inserted  into  the 
nozzle  assembly  through  a  hole  bored  in  the  vaned  elbow  (see 
Figure  2.4)  and  positioned  on  the  center  line  of  the  flow.  The  tip 
of  the  electrode  was  adjusted  to  be  approximately  4  nozzle  diameters 
upstream  of  the  nozzle  exit.  A  length  of  insulated  wire  connected  the 
electrode  to  the  negative  terminal  of  a  D.C.  power  supply.  The  other 
terminal  was  grounded  to  the  tunnel.  In  this  configuration,  hydrogen 
bubbles  were  generated  at  the  electrode.  By  varying  the  current  to 
the  electrode,  different  nuclei  populations  in  the  flow  were  artifi¬ 
cially  created. 
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2.6  Inception  Detection  System 

Figure  2.6  is  a  diagrammatic  layout  of  the  detection  system  em¬ 
ployed  in  the  present  studies.  The  sensor  was  a  microphone  pressed 
against  the  window  of  the  tunnel  and  placed  directly  above  the  cavita- 
ting  region.  The  microphone  was  insulated  from  surrounding  noises. 

On  cavitation  inception,  the  noise  generated  by  the  exploding  and  col¬ 
lapsing  cavities  was  picked  up  by  the  microphone.  The  signal  from  the 
microphone  was  then  fed  through  the  preamplifier  section  of  a  com¬ 
mercial  tape  deck.  From  the  preamplifier,  the  amplified  signal  was 
band  passed  by  a "KRONHITE" band  pass  filter  set  in  the  range  of  1.9  to 
2.1  kHz.  The  output  from  the  tape  was  also  monitored  by  a  pair  of 
headphones  to  ensure  that  the  signal  picked  up  by  the  microphone  was 
that  of  a  cavitating  jet  at  inception  and  not  of  some  other  sources. 

The  filtered  signal  was  further  amplified  by  a  variable  gain  amplifier. 
In  the  present  work,  the  gain  of  the  amplifier  was  set  at  7.  The 
signal  was  finally  fed  into  a  frequency  counter.  When  a  count  of  5 
per  sec  or  more  was  noted,  inception  was  called.  All  inception  data 
were  obtained  this  way. 

In  addition  to  measuring  the  incipient  cavitation  number  under 
different  flow  conditions,  it  was  desired  to  observe  the  flow.  Flow 
visualization  was  accomplished  with  the  aid  of  the  holocamera  and  the 
schlieren  system.  This  equipment  will  now  be  described. 

2.7  Holocamera  and  Reconstruction  System 

The  present  holographic  system  is  a  completely  new  design  by 
Katz  The  various  components  in  the  holocamera  are  shown  in 
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Figure  2.1.  The  principal  unit  is  a  ruby  rod  measuring  7.62  cm 
(3.00  inch)  in  length  and  0.6  cm  (0.25  inch)  in  diameter.  The  ends 
of  the  ruby  rod  are  highly  polished  plane  surfaces.  Excitation  of 
the  ruby  is  brought  about  by  a  helical  xenon  flashlamp.  Cylindrical 
aluminum  reflectors  which  cover  the  flashlamp  are  used  to  reflect  the 
light  from  the  flashlamp  into  the  ruby  rod.  Cooling  of  the  ruby  rod- 
flashlamp  is  accomplished  by  continuously  recirculating  deionized  water 
through  this  assembly  during  operation. 

The  optical  cavity  is  formed  by  placing  two  flat  surface  mirrors  on 
both  sides  of  the  ruby  rod-xenon  flashlamp  assembly.  They  are  held  paral¬ 
lel  to  each  other.  The  back  mirror  (component  3  in  Fig. (2.1))  is  a  100$ 
reflecting  dielectric  surface,  while  the  front  mirror  is 
sapphire  ethalen  and  is  only  60$  reflecting.  The  front  mirror  also 
assists  in  the  longitudinal  mode  control  of  the  laser.  Transverse 
mode  control  is  achieved  by  two  iris  apertures  (components  4  and  13) 
placed  in  the  optical  cavity. 

When  the  flashlamp  is  triggered,  about  1  kJ  of  energy  is  dis¬ 
charged  through  it.  When  the  optical  cavity  is  open,  the  resulting 
light  emitted  by  the  ruby  rod  races  back  and  forth  between  the  two 
mirrors  and  at  each  passing,  more  electrons  in  the  ruby  are  excited. 

Hence  a  population  inversion  is  created  and  eventually  results  in  the 
lasing  of  the  ruby  rod.  Since  the  duration  that  the  flashlamp  is 
on  is  much  longer  than  that  for  a  single  lasing  of  the  ruby,  numerous 
light  pulses  are  generated  in  a  single  flash  of  the  flashlamp.  This 
is  quite  undesirable  because  it  causes  multiple  exposure  to  the 
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holographic  film.  For  this  reason,  the  cavity  has  to  be  "Q  switched". 

★ 

"Q  switching"  is  accomplished  by  a  KD  P  cylindrical  ring  electrode 
pockel's  cell  and  a  calcite  glan  polarizer.  This  process  works  as 
follows: 

Light  emitted  by  the  ruby  is  first  polarized  by  the  glan  polarizer. 
★ 

It  then  enters  the  KD  P  crystal.  A  high  voltage  is  applied  to  the 
crystal  so  that  the  angle  of  polarization  of  the  incoming  light  is 
turned  through  45°.  After  exiting  from  the  pockel's  cell,  the  light 
is  reflected  back  by  the  mirror  into  the  cell  again.  On  this  second 
passage  through  the  cell,  the  polarization  angle  is  rotated  through 
another  45°.  When  the  reflected  light  reaches  the  polarizer,  it  is 
now  rotated  through  90°  to  its  original  polarization  angle  and  thus 
not  permitted  into  the  ruby.  Population  inversion  is  thus  avoided  and 
no  lasing  occurs.  In  short,  when  the  pockel’s  cell  is  on,  the  optical 
cavity  is  closed. 

When  the  voltage  to  the  pockel  cell  is  dropped  to  zero,  the 
cavity  is  effectively  open  to  light  oscillations  which  results  in 
the  eventual  lasing  of  the  ruby.  Hence  the  idea  is  to  drop  the  voltage 
to  the  cell  for  a  sufficiently  long  time  during  the  flashing  of  the 
xenon  lamp  so  that  the  ruby  is  allowed  to  lase  only  once  (or  twice,  if 
double  pulsing  is  desired).  This  is  accomplished  by  setting  a  delay 
time  on  the  power  source  to  the  pockel’s  cell.  The  delay  time  here  is 
defined  as  the  time  that  elapsed  between  the  instant  the  flashlamp  is 
triggered  to  the  instant  the  voltage  to  the  cell  is  dropped  to  zero. 

The  right  setting  of  the  delay  time  is  obtained  through  trial  and 
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error  and  Is  usually  around  800  us.  It  should  be  pointed  out  that  the 
duration  of  the  laser  pulse  is  between  20  to  50  nanoseconds. 

A  beam  splitter  directs  about  4%  of  the  beam  emitted  from  the 
cavity  into  a  pin  diode  which  monitors  the  output  of  the  laser. 

The  remaining  beam  passes  through  a  microscope  objective  which  expands 
the  beam.  The  resulting  beam  is  spatially  filtered  by  an  11  jm  pin 
hole  and  finally  collimated.  The  collimated  beam  is  used  to 
illuminate  the  sample  volume.  A  camera  placed  on  the  opposite  window 
of  the  tunnel  is  used  for  recording  the  images.  The  recording  film 
used  is  70  mm  Agfa  Gevaert  film. 

Holograms  obtained  are  then  reconstructed  in  the  reconstruction 
system.  This  system  is  illustrated  in  Figure  2.7.  The  hologram  to 
be  examined  is  mounted  on  an  x-y-z  vernier  carriage  and  is  illuminated 
by  a  5  mW  Spectra  Physics  He-Ne  laser.  The  beam  is  expanded  and 
collimated  in  the  same  manner  as  in  the  holocamera.  The  reconstructed 
image  is  magnified  220  times  with  a  microscope  objective  and  imaged  on 
a  videcon.  This  magnified  image  is  then  displayed  on  a  closed 
circuit  TV  monitor.  The  sizes  of  the  nuclei  are  read  directly  off  the 
screen  of  the  monitor  with  the  aid  of  a  reticle. 

2.8  Schlleren  System 

A  comprehensive  description  of  the  present  schlieren  system  is 
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found  in  the  thesis  by  Gates  v  A  schematic  drawing  of  this  system 


is  shown  in  Figure  2.8.  Basically,  the  light  source  was  generated  by 
discharging  two  0.05  microfarad,  10  kV  capacitors  through  a  tungsten 
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electrode.  The  duration  of  the  spark  was  about  two  microseconds  thus 
permitting  photographs  of  the  instantaneous  flow  field  under  observa¬ 
tion  to  be  recorded.  The  light  emitted  from  the  spark  source  was 
collimated  by  a  12  inch  focal  length  lens  and  was  directed  across  the 
test  section  of  the  water  tunnel.  In  the  present  work,  the  required 
"schlieren  effect"  was  created  by  heating  the  water  in  the  jet  to 
about  2°  to  4°  C  above  the  temperature  of  the  surrounding  fluid.  It 
is  worth  mentioning  that  preliminary  work  showed  that  a  jet  heated 
to  only  1°  C  above  the  ambient  liquid  could  also  be  photographed. 
However,  poor  contrast  on  the  recording  film  was  observed  especially 
at  distances  far  downstream  of  the  jet  exit.  Photographs  of  "Epsom 
salt  solution"  jets  (0.5  percent  concentration)  were  also  taken  when 
high  contrast  was  desired. 

The  scattered  light  emerging  from  the  test  section  was  then 
focused  onto  a  knife  edge.  The  edge  was  adjusted  to  block  a  major 
portion  of  the  beam.  Two  knife  edge  positions  were  used:  (i)  with 
the  cutting  edge  parallel  to  the  direction  of  flow  and  (ii)  with  the 
edge  perpendicular  to  the  flow  direction.  The  image  was  finally 
recorded  on  Kodak  4x5 "Royal -X-Pan" film.  Due  to  the  size  limitation 
of  the  recording  film,  most  of  the  photographs  were  taken  at  low 
image  magnification  (approximately  3)  to  enable  a  larger  portion  of 
the  jet  to  be  recorded. 
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Fig.  2.2(b).  Photograph  of  the  dearation  tank. 


Fig.  2.3.  Graph  displaying  the  performance  curve  of  the  deaeration  system. 
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Fig.  2.4.  Diagram  showing  the  geometrical  details  of  the  nozzle  assembly. 


Fig.  2.7.  Diagram  of  the  holographic  reconstruction  system. 
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3.  EXPERIMENTAL  PROCEDURES  AND  MEASUREMENT  TECHNIQUES 

3.1  Introduction 

The  experimental  procedures  followed  during  the  course  of  the 
present  studies  are  described  in  this  chapter.  The  experimental  ob¬ 
jectives  called  for  the  measurement  of  the  incipient  cavitation 
indices  under  different  flow  conditions.  For  each  nozzle  size,  the 
cavitation  index  was  measured  as  a  function  of  (a)  the  dissolved  air 
content  of  the  water,  (b)  the  exit  velocity  of  the  jet  and  (c)  the 
nuclei  population  in  the  jet.  The  nuclei  population  in  the  jet 
was  altered  by  seeding  the  flow.  Three  nozzle  sizes  were  investigated. 
They  were  3.17  run,  4.76  nm  and  6.35  inn  in  diameter. 

Flow  visualization  was  also  carried  out.  This  was  accomplished 
by  holography  and  schlieren  photography.  Furthermore,  the  flow  was 
observed  under  strobe  light.  In  addition,  the  effect  of  the  dis¬ 
solved  air  content  on  the  nuclei  population  density  was  studied  with 
the  aid  of  the  holocamera.  Prior  to  the  actual  data  acquisition  and 
flow  observation,  several  preliminary  preparations  had  to  be  performed. 
These  preliminary  steps  taken  are  described  first,  followed  by  a 
description  of  the  measurement  techniques  employed. 

3.2  Preliminary  Preparations 

The  nozzle  assembly  was  mounted  on  the  top  window  of  the  test 
section  of  the  LTWT  as  shown  in  Figure  3.1.  It  was  positioned  such 
that  the  axis  of  the  nozzle  was  on  the  center  line  of  the  water  tunnel. 
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The  tungsten  electrode  used  for  "electrolysis  seeding"  was  inserted 
into  the  elbow  of  the  nozzle  assembly,  and  was  adjusted  until  the  tip 
was  about  four  nozzle  diameters  from  the  lip  of  the  nozzle.  The 
tunnel  water  was  then  filtered  for  six  hours.  This  was  followed  by 
deaerating  the  water  to  the  desired  test  condition.  The  microphone 
used  for  cavitation  detection  was  placed  against  the  window  of  the 
tunnel  such  that  it  was  directly  above  the  cavitating  regime  of  the 
jet. 

On  completion  of  these  preliminary  preparations,  the  data  ac¬ 
quisition  process  began. 

3.3  Measurement  of  the  Incipient  Cavitation  Index 

Some  water  from  the  water  tunnel  was  pumped  into  the  pressure 
vessel.  Pressurization  of  the  vessel  was  accomplished  with  the 
Institute  compressed  v  ir  supply.  The  inlet  pressure  to  this  vessel 
was  regulated  witn  a  pressure  regulator.  A  jet  of  water  was  then  dis¬ 
charged  into  the  stationary  tunnel  sitting  at  an  ambient  pressure  of 
approximately  0.9  atmospheres  (13  psia).  The  pressure  in  the  tunnel 
was  then  quickly  lowered  until  cavitation  inception  in  the  jet 
occurred.  As  explained  earlier,  inception  was  called  at  the  instant 
5  events  per  second  or  more  were  indicated  by  the  detection  system. 
The  pressure  in  the  tunnel  and  the  exit  velocity  of  the  jet  at 
inception  were  noted.  From  these  measurements,  the  incipient  cavita¬ 
tion  index  was  calculated  from  the  expression  below: 
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where  P^  =  the  pressure  in  the  tunnel  when  inception  was  called 

P y  =  vapor  pressure  of  water  evaluated  at  the  bulk  temperature 
of  the  water 
p  =  density  of  water 
and  V.  *  exit  velocity  of  the  jet. 
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The  tunnel  pressure  was  raised  to  its  initial  value  and  allowed 
to  sit  at  that  pressure  for  about  30  seconds  before  similar  runs  were 
performed. 

The  incipient  index  was  also  measured  when  the  flow  was  seeded. 
Seeding  was  done  by  electrolysis  of  the  water  and  the  following 
currents  were  used  to  create  different  nuclei  populations:  0.05  mA, 

0.1  mA,  0.2  mA,  0.4  mA,  0.6  mA,  1.0  mA,  2.0  mA,  and  5.0  mA. 

All  the  above  tests  were  repeated  over  a  range  of  jet  velocities. 
Different  velocities  in  the  jet  were  realized  through  pressurization  of 
the  vessel  to  different  values  and  not  by  flow  throttling. 

For  each  nozzle  size,  the  above-mentioned  tests  were  carried  out 
sequentially  with  the  tunnel  water  at  different  levels  of  air  satura¬ 
tion.  Jet  water  with  the  following  dissolved  air  content  was  investi¬ 
gated:  14.1  ppm,  10.9  ppm,  7.7  ppm  and  4.2  ppm.  (Saturation  is  at 
15  ppm.) 


The  incipient  cavitation  index  corresponding  to  a  particular  flow 
condition  was  obtained  by  averaging  the  results  of  ten  similar  runs. 

3.4  Flow  Visualization  and  Nuclei  Measurements 

The  flow  was  observed  under  stroboscopic  light  as  the  jet  was 
brought  to  cavitation  inception  and  beyond.  In  addition,  schlieren 
photography  was  also  employed  to  give  a  permanent  record  of  the 
flow  observation.  Two  sets  of  schlieren  photographs  were  obtained. 

The  first  set  was  taken  with  the  cutting  edge  of  the  cutoff  plate 
parallel  to  the  direction  of  flow,  while  the  second  set  was  with  the 
edge  perpendicular.  Although  these  photographic  records  gave  valuable 
information  such  as  the  spread  rate  of  the  jet  and  the  location  of 
the  cavities  of  inception,  they  provided  no  clue  to  the  miroscopic 
nuclei  population  in  the  jet.  This  is  a  major  setback  since  such 
information  would  definitely  give  a  better  insight  into  the  cavita¬ 
tion  phenomena  in  jets  on  a  whole.  This  shortcoming  was  overcome 
with  holographic  studies. 

As  a  result  of  some  preliminary  work,  it  was  found  that  with  the 
present  in-line  holography  setup,  holograms  of  heated  jets  could  not 
be  reconstructed  and  studied.  This  was  because  the  reference  beam 
itself  was  scattered  by  the  density  field,  even  though  the  water  in 
the  jet  was  heated  to  only  1/2°  C  above  the  surrounding  water.  (Ooi 
and  Katz  Consequently  only  holograms  of  unheated  jets  were 


obtained  and  analyzed. 
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The  three  flow  visualization  techniques  employed  above  were  per¬ 
formed  for  all  the  flow  conditions  previously  mentioned  in  Section  3.3. 

A  new  method  was  devised  to  measure  the  pressure  fluctuation  in 
the  present  jets.  The  motivation  for  this  study  and  the  measurement 
techniques  used  are  presented  in  the  next  chapter. 


Photograph  showing  the  counted  nozzle  assembly  on  the  top  window  of  the  water  tunnel. 
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4.  MEASUREMENT  OF  PRESSURE  FLUCTUATIONS  IN  WATER  JETS 
4.1  Introduction 

Since  the  pioneering  work  of  Rouse  the  study  of  cavitation 

in  water  jets  have  been  actively  pursued.  The  cavitation  inception 
nimbers  have  been  measured  by  numerous  investigators,  for  example, 

Hoyt  ^8)  t  Lienhard  and  Goss  and  Jorgensen  ^63^.  Although  the 
measured  indices  were  quantitatively  different,  the  results  of  these 
researchers  nevertheless  indicated  that  cavitation  inception  always 
occurred  at  a  mean  pressure  which  was  considerably  higher  than  the 
vapor  pressure  of  the  liquid.  Obviously,  this  observed  fact  was  in 
direct  contradiction  to  the  long  held  belief  that  cavitation  would 
only  occur  if  the  pressure  is  at  or  below  the  vapor  pressure.  To 
explain  this  discrepancy,  it  was  postulated  that  local  pressures  at 
numerous  positions  in  the  jet  dip  to  or  below  the  vapor  pressure  due 
to  pressure  fluctuations.  This  will  consequently  give  rise  to  cavita¬ 
tion  inception.  If  this  is  indeed  the  case,  the  departure  of  the 
instantaneous  local  pressure  from  the  mean  would  be  enormous. 

One  of  the  first  attempts  to  experimentally  determine  the  pres¬ 
sure  fluctuations  in  axisymmetric  jets  was  carried  out  by  Rouse 
His  work  was  later  extended  by  Sami  et_al_  Barefoot 

made  similar  measurement  and  although  his  measured  values  were  sub¬ 
stantially  different  from  those  reported  by  Sami  et  al_  the 

data  followed  the  same  general  trend.  Other  workers  who  have  reported 
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work  on  pressure  fluctuations  or  spectra  measurements  include  Fuchsv  ’ 
Arndt  and  George  and  Lau  et^a^  ^6).  the  above 
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measurements  were  made  In  atr  jets  by  introducing  a  pressure  trans¬ 
ducer  or  a  mtcropfione  Into  the  flow\ 

The  highest  rms  pressure  fluctuations  ever  recorded  was  about 
22  percent  of  the  dynamic  head  of  the  jet  (Sami  et  al  These 

measurements  are  themselves  Insufficient  to  account  for  cavitation 
inception  in  some  jets,  where  It  has  Been  reported  that  the  measured 
incipient  index  could  Be  as  high  as  3  in  a  13  cm  jet. (Lienhard  and 
Goss  In  other  jets,  the  occurrence  of  inception  could  only 

be  explained  if  the  peak  pressure  fluctuations  are  at  least  several 
factors  higher  than  the  measured  rms  values.  Arndt  and  George 
have  stated  that  the  peak  to  peak  values  of  pressure  fluctuations 
could  Be  as  much  as  ten  times  the  measured  rms. 

In  view  of  the  above  discussion,  it  is  understandable  that 
there  is  a  need  to  determine  if  such  high  pressure  fluctuations  occur 
in  the  jets  that  are  presently  investigated.  Since  these  jets  are 
relatively  small  in  size  (the  largest  is  only  6.35  nin},  the  normal 
approach  of  measuring  with  a  pressure  transducer  is  not  suitable. 

Hence  there  is  a  need  to  devise  a  way  to  measure  the  pressure  fluc- 

tuati ons . 

It  is  known  that  when  a  bubble  is  subjected  to  a  pressure  field, 
it  will  deform  in  shape  and  volume  accordingly.  If  this  bubble  is 
designed  to  meet  certain  requirements,  which  will  be  considered 
shortly,  it  can  serve  as  a  good  pressure  sensor.  In  principle,  the 

pressure  field  in  the  jet  could  then  be  mapped  out  by  following  the 

bubble  as  it  traverses  the  jet.  The  theory  behind  the  idea  of  using 
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such,  a  buhhle  as  a  static  pressure  sensor  ts  considered  In  the  next 
section. 

4,2  Theory  Behind  Stattc  Pressure  Measurements  Using  Bubbles  as 
Pressure  Sensors 

The  equation  governing  the  motfon  of  a  spherical  Bubble  is  given 
by  the  familiar  Rayletgh-Plesset  equation  (PTesset  and  Prosperetti 
viz. 
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Rewriting,  we  have 


t 

pjti  =  [pgcTj  +  pvcti  -  ir  nr]  -  p[R* +  3/2^2]  o ) 

For  a  given  bubble  size,  the  static  pressure  in  the  jet  at  any  instant 

•  »  « 

can  then  be  computed  from  the  above  equation  once  R  and  R  are  known. 
In  addition,  if  the  bubble  is  designed  to  have  a  response  time  that 
is  short  in  comparison  to  the  time  scale  of  the  typical  pressure 
fluctuations, then  the  inertia  and  the  viscous  terms  in  the  Rayleigh- 
Plesset  equation  can  be  neglected  as  shown  in  Appendix  1.  Equation 
(1)  then  simplifies  to; 

Pjt)  *  PgCT)  +  PV(T)  -  .  12) 
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In  theory,  the  Instantaneous  pressure  at  various  positions  in  the 
jet  can  then  be  computed  from  a  continuous  photographic  record  of  the 
size  of  a  bubble  as  it  crosses  the  jet  from  a  reference  point.  In 
practice  this  is  more  difficult  to  achieve  since  the  events  are 
occurring  at  a  high  speed,  Furthermore.even  if  somehow  one  succeeded 
in  tracing  the  path  of  the  bubble  under  observation,  the  pressures 
indicated  by  the  bubble  correspond'  to  different  instants  in  time. 
Consequently,  one  would  not  have  the  true  view  of  the  pressure  field 
in  the  jet  at  any  one  instant. 

Consider  a  string  of  bubbles  always  of  the  same  size  at  the 

point  of  injection,  which  is  being  continuously  introduced  into  the 

jet  at  a  reference  position.  An  instantaneous  photograph  would  then 

reveal  the  sizes  of  the  bubbles  at  various  positions  in  the  jet.  Any 

change  in  the  bubble  size  from  the  reference  size  could  be  interpreted 

to  be  due  to  deviation  of  pressures  from  the  mean.  I -f  P  is  taken 

9o 

to  be  pressure  of  the  gas  in  the  bubble  at  the  point  of  injection, 
then  the  pressure  of  the  gas  at  some  position,  x,  is  given  by 


3y 


where  R„  and  R  are  the  radii  of  the  bubble  at  the  reference  point 
ox 

and  at  position  x  ,  respectively,  and  y  is  the  polytropic  index. 
The  instantaneous  pressure  at  point  x  can  then  be  computed  from: 


(3) 


In  tMs  manner,  an  instantaneous  record  of  the  presssure  field  across 
the  entire  jet  could  6e  determined.  Furthermore,  this  measuring 
technique  has  the  added  advantage  that  the  "probes1'  do  not  cause  a 
significant  flow  disturbance. 

As  mentioned  earlier,  the  bubbles  can  serve  as  good  pressure 
sensors  only  if  they  are  tailored  to  meet  certain  requirements.  These 
requirements  are  summarized  below. 

4.3  Summary  of  Requirements 

It  has  been  mentioned  that  the  technique  adopted  required  a 
string  of  bubbles  of  a  uniform  referenced  size.  Also,  in  order  for 
the  bubbles  to  be  good  indicators  of  the  instantaneous  pressure,  there 
should  be  no  lag  time  in  the  response  to  a  sudden  pressure  change. 

This  implies  that  the  response  time  of  the  bubble  should  be  at  least 
an  order  of  magnitude  faster  than  the  typical  time  scale  of  the  pressure 
fluctuations. 

The  bubbles  are  also  required  to  be  small  enough  so  that  the 

pressures  they  see  are  a  good  approximation  of  the  local  "point" 

pressures.  However,  they  should  not  be  so  small  that  the  change  in 

(431 

volume  could  not  be  easily  and  accurately  determined.  Fuchs  ^  , 

Michalke  and  Fuchs  and  Lau  e£  al_  were  among  the  relatively 

few  who  have  measured  the  pressure  spectra  in  jets.  Their  work 
demonstrated  that  the  dominant  portion  of  the  spectra  was  shifted  to 
lower  frequencies  with  increasing  distance  from  the  nozzle  exit. 

The  results  of  the  three  groups  of  investigators  mentioned  above  also 
showed  that  the  normalized  power  spectra  peaked  at  a  Strouhal  number 
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(defined  as  1  of  0.45,  and  beyond  that  the  spectra  dropped  rapidly 
to  insignificant  values,  ffarefoot  also  reported  comparable 
results.  This  meant  that  in  the  present  jets,  the  frequency  of  the 
dominant  portion  of  the  pressure  fluctuations  is  of  the  order  of  a 
kilohertz  which  translated  into  a  time  scale  of  one  millisecond. 

Arndt  and  George  have  derived  an  equation  for  a  character¬ 
istic  bubble  response  time  namely 


pr; 
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where  RQ  is  the  radius  of  the  bubble  and  S,  the  surface  tension. 

Therefore  for  a  bubble  of  100  ym  radius,  the  response  time  is  of  the 
order  of  a  tenth  of  a  millisecond.  This  is  an  order  of  magnitude  shorter 
than  the  time  scale  for  the  pressure  fluctuations.  In  addition,  any 
change  in  volume  of  the  bubble  caused  by  the  fluctuation  of  pressures 
can  be  easily  detected.  Hence,  a  100  ym  bubble  was  taken  to  be  an  ideal 
size  for  the  present  investigation. 

With  the  preferred  bubble  size  determined,  steps  were  taken  to 
create  these  bubbles.  This  is  discussed  in  the  next  section. 

4.4  Preliminary  Work 

It  was  believed  that  the  desired  bubble  size  could  be  generated 
by  electrolysis.  Trials  carried  out  In  a  beaker  of  water  using  an 
electrode  which  was  similar  in  description  to  that  employed  in  the 
flow  seeding  process  have  shown  that  the  size  of  the  bubbles  was 
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determined  fay  the  magnitude  of  the  electrolysis  current,  For  example 
when  the  current  was  5  mA,  faufafales  as  Big  as  1  mm  in  size  were  observed. 
Consequently  we  thought  that  the  required  bubble  size  could  be  achieved 
by  a  suitable  combination  of  the  flow  velocity  and  the  electrolysis 
current.  This  desired  condition,  of  course,  was  to  Be  obtained 
through  a  trial  and  error  process. 

However,  when  the  tests  were  carried  out  in  the  nozzle  with  the 
electrode  inserted  in  the  elbow  as  before,  it  was  found  that  the 
size  of  the  generated  bubbles  when  there  was  a  flow  past  the  electrode 
was  only  about  10  vm.  Furthermore,  neither  the  magnitude  of  the  flow 
velocity  nor  the  electrolysis  current  appeared  to  have  any  major 
effect  on  the  size  of  these  electrolyzed  hydrogen  bubbles.  Different 
sizes  and  shapes  of  the  tips  of  the  electrodes  were  tried  too  but  to 
no  avail.  It  was  suggested  that  a  larger  size  of  bubbles  could  be 
obtained  if  the  electrode  was  placed  normal  to  the  direction  of  flow. 
However,  based  on  Kuiper's  work  it  was  inferred  that  a  relatively 

broad  spectrum  of  bubble  sizes  would  be  produced  in  such  a  configura¬ 
tion.  This  is  undesirable  since  we  want  only  one  particular  referenced 
size.  It  was  finally  decided  to  try  a  new  approach. 

Air  injection  was  considered  next.  A  length  of  stainless  steel 
tubing  with  an  internal  diameter  of  100  ym  was  used  as  the  injector. 

It  was  found  that  although  the  bubbles  created  were  approximately  of 
the  right  size,  they  were  not  very  spherical.  This  was  because  the 
orifice  at  the  exit  of  the  tubing  was  not  perfectly  circular,  as  a 


is  • 


result  of  it  being  partially  squashed  during  the  cutting  of  the  tubing 
to  the  prescribed  length.  Since  thfs  shortcoming  could  not  be  easily 
avoided,  it  was  decided  that  a  circular  glass  nozzle  be  used  as  the 
injector  instead.  This  nozzle  was  made  By  drawing  a  small  glass 
tubing  down  to  the  necessary  dimensions.  This  method  eventually 
proved  to  be  successful.  The  bubble  generator  and  the  accompanying 
injection  system  is  described  below. 


4.5  Description  of  the  Bubble  Generator  and  the  Injection  System 

A  diagram  of  the  injection  system  is  illustrated  in  Figure  4.1, 
Compressed  air  supply,  regulated  to  4  psig,  was  filtered  before 
being  fed  to  the  bubble  generator  to  prevent  clogging  of  the  small 
orifice  of  the  glass  nozzle.  A'^uprcf* meter  valve  C"S  series")  located 
downstream  of  the  air  filter  allowed  the  small  air  flow  rate  through 
the  system  to  be  adjusted  until  a  tiny  stream  of  bubbles  was  seen. 

A  shutoff  valve  was  also  provided  to  stop  the  flow  of  air  into  the 
bubble  generator  in  between  runs,  without  having  to  tamper  with  the 
setting  of  the  meter  valve,  A  by-pass  line  permitted  any  back  flow 
of  liquid  into  the  bubble  generator  to  be  quickly  expelled  prior  to 
an  actual  run. 

The  bubble  generator  consisted  of  a  glass  nozzle  drawn  from  a 

1.6  mm  O.D.  glass  tubing.  The  inner  diameter  at  exit  was  designed 
to  be  approximately  100  pm  CO .004"},  The  glass  nozzle  was  attached 
to  a  length  of  brass  tubing  and  was  held  together  by  "intra  medic 
tubing"  as  shown  in  Figure  4,2,  This  whole  unit  was  introduced  into 
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the  jet  In  the  same  manner  as  for  the  electrode  used  in  electrolysis 
bubble  generation.  This  has  been  described  in  Section  2,4  of  Chapter 
2, 

The  tip  of  the  bubble  generator  was  positioned  on  the  center 
line  at  approximately  four  diameters  upstream  of  the  nozzle  exit.  A 
length  of  1,6  nin  I.D,  "intra  medic  tubing"  connected  the  other  end 
of  the  bubble  generated  to  the  air  supply.  This  completed  the  air 
injection  system. 

4.6  Test  Procedures 

For  a  given  jet  velocity,  the  air  flow  rate  through  the  injector 
was  adjusted  until  a  stream  of  individual  bubbles  was  seen  under 
stroboscopic  lighting.  A  hologram  of  the  flow  field  was  then  taken. 

In  this  way,  a  series  of  holograms  were  obtained.  The  procedure  was 
repeated  for  all  three  nozzle  sizes.  It  should  be  mentioned  at  this 
point  that  the  bubbles  generated  were  actually  less  than  100  ym,  but 
we  shall  still  assume  that  they  are  of  this  size  in  our  argument. 

In  all  the  above  tests,  the  ambient  pressure  was  kept  at  approxi¬ 
mately  16  psia.  The  tests  were  carried  out  at  a  jet  exit  velocity  of 
13  m/s  (=  40  ft/s} .  This  velocity  was  chosen  because  the  resulting 
pressure  fluctuations  would  be  large  enough  to  be  easily  detected. 

On  the  other  hand,  the  fluctuation  would  not  be  so  intense  to  cause 
cavitation  In  the  jet.  All  the  measurements  were  made  using  water 
with  a  dissolved  air  content  of  approximately  11  ppm  (saturation  Is 
at  15  ppm}. 


The  sizes  of  the  buhhles  as  a  function  of  the  jet  position  was 
then  measured  from  the  reconstructed  images  of  the  holograms  and 
tabulated. 

4,7  Data  Reduction 

Before  the  instantaneous  pressures  at  various  locations  in 
the  jet  could  be  computed  from  the  tabulation  of  bubble  sizes,  two 
aspects  of  the  flow  had  to  be  considered.  The  first  of  these  had 
to  do  with  the  “ringing"  effect  of  the  bubbles  as  they  were  accelerated 
through  the  flow  nozzle.  The  other  dealt  with  the  rate  at  which  the 
injected  bubbles  were  dissolving  into  the  water.  Both  of  these  will 
be  considered  in  turn. 

Following  an  argument  similar  to  that  presented  by  Franklin  and 
McMillan  it  can  be  assumed  that  the  motion  of  the  bubble  after 

exiting  from  the  flow  nozzle  is  governed  by  the  equation: 

RCtl  *  Re  e"6t  cos  t 

where  Rg  =  initial  radius  of  the  bubble  after  exiting  from  the  nozzle, 

6  *  damping  constant 

u)Q  *  natural  frequency  of  the  bubble  and 

T  a  V  ’ 

Devin  has  determined  theoretically  that  the  damping  con¬ 
stant  for  a  100  ym  bubble  Is  0.1.  He  showed  that  his  results  compare 
favorably  with  the  experimental  findings  of  several  investigators 
although  his  theoretical  value  was  slightly  lower.  The  same  value 
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for  g  was  also  reported  by  Clay  and  Medwtn  As  a  conservative 

estimate,  the  numerical  value  of  g  was  thus  taken  to  be  Q,1  In  the 

present  work.  Hence  the  time  that  elapsed  before  the  injected  bubbles 

-fit 

acquired  steady  state  could  be  determined  from  a  plot  of  e  p  cos  t 
vs  t.  Such  a  plot  Is  shown  in  Figure  4.3.  It  can  be  inferred  from 
this  plot  that  the  oscillations  would  have  died  down  for  t  >_  40, 

This  implied  that  equilibrium  was  reached  after  about  30  us.  That 
is,  at  a  jet  exit  velocity  of  13  m/s,  the  oscillations  would  have 
virtually  stopped  at  distances  greater  than  about  0,5  mm.  This  has 
been  confirmed  from  holographic  records  of  the  injected  bubbles  near 
the  nozzle  lip  which  showed  no  evidence  of  any  oscillations  in  the 
string  of  bubbles.  A  photographic  evidence  of  this  is  shown  in 
Figure  4.4. 

We  shall  now  direct  our  attention  to  the  question  of  the  rate 
at  which  the  injected  bubbles  were  dissolving.  We  would  expect  a 
huge  error  in  the  measured  volumes  of  the  bubbles  if  the  bubbles 
were  going  into  solution  at  a  rate  comparable  to  the  time  they  took 
to  traverse  the  region  under  observation.  By  interpolating  the 
results  of  Epstein  and  Plesset  ^4)  we  have  found  that  it  would  take 
about  460  seconds  for  a  100  um  bubble  to  completely  dissolve  in  water 
that  Is  75%  saturated  (a  dissolved  air  content  of  11  ppm).  In  addi¬ 
tion,  it  has  been  estimated  that  the  traverse  time  of  the  bubble 
in  the  current  Investigation  was  of  the  order  of  milliseconds.  Hence 
it  can  be  concluded  that  the  dissolved  volume  of  the  bubbles  would  be 
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negligible  and  can  be  ignored, 

As  was  discussed  in  Section  4,2,  the  Instantaneous  pressure  at 
any  position  is  given  By 

\Ctl  ’  *  PV  '  f 

Since  dry  air  was  used  in  the  present  tests,  the  vapor  pressure  inside 
the  bubble  was  neglected,  He  believe  this  is  justified  because  in 
the  extremely  short  traverse  time  involved,  little  or  no  evaporation 
was  expected  to  occur.  Thus  equation  03 1  reduces  to 
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It  can  be  seen  from  equation  051  above  that  once  P  is  known  then 

9o 

P^  Ct)  at  any  position  could  be  calculated  using  the  tabulation  of  the 
x 

bubble  sizes. 

The  injected  bubbles  have  been  shown  to  acquire  steady  state 

conditions  very  soon  after  emerging  from  the  flow  nozzle.  Therefore, 

the  initial  gas  pressure  P  can  be  obtained  from  the  reduced  Rayleigh- 

9o 

Plesset  equation,  that  is 


P  m  p  +  jpL 

9n  0  K 
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where  PQ  »  mean  ambient  pressure  In  the  jet  and 

Rq  »  the  radius  of  the  Injected  BuBBle  at  the  reference  position. 

On  substitution  of  P  into  Equation  C3i»  the  equation  becomes 

9o 

p”xctl '  cp. * 11  ^  f  t7) 

Observations  of  reconstructed  images  of  holograms  revealed  that 
in  the  region  beyond  the  potential  core  of  the  jet,  some  of  the  bub¬ 
bles  did  not  retain  their  spherical  forms  but  were  ellipsoidal  in 
shape.  An  example  of  this  situation  is  shown  in  Figure  4.5,  The 
question  that  came  to  mind  at  this  point  was  whether  these  bubbles 
still  served  as  good  pressure  sensors.  To  answer  this  question,  we 
have  to  rely  on  the  findings  of  Strasberg  ^08^,  He  reported  that 
the  natural  frequency  of  an  ellipsoidal  bubble  could  be  very  accurately 
estimated  from  the  response  time  of  a  spherical  bubble  of  the  same 
volume.  His  calculations  showed  that  the  natural  frequency  of  bubbles  of 
this  shape  was  only  slightly  dependent  on  the  ratio  of  the  major  to 
the  minor  axis  of  the  ellipsoid.  In  the  present  measurements,  the 
largest  change  in  volume  ever  detected  was  about  3.5  times  the 
original.  This  translated  Into  an  increase  of  about  50  percent  in  the 
Initial  bubble  size  or  a  bubble  of  150  ym  radius.  It  can  easily  be 
shown  from  Equation  (4)  that  a  bubble  of  this  size  still  has  a 
relatively  short  response  time.  Hence  it  can  be  concluded  that  all 
the  bubbles  remained  as  good  Indicators  of  the  Instantaneous  pressure. 
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In  order  to  determine  the  volume  of  these  ellipsoidal  bubbles, 
both  the  semi -mi nor  and  semi -major  axis  need  to  be  measured.  However, 
it  is  not  possible  to  get  the  third  dimension  of  the  bubble  from  the 
reconstructed  image.  It  was  thought  that  this  dimension  could  be 
determined  theoretically.  However,  it  has  been  pointed  out  by 
Barthds-Biesel  and  Acrivos  that  "the  problem  of  determining  theore¬ 
tically  the  shape  of  a  single  droplet  freely  suspended  in  an  unbounded 
incompressible  liquid  undergoing  a  shearing  motion  is  very  complex  and 
no  general  solution  is  at  present  available."  Most  of  the  work  done 
in  this  area  has  been  on  very  simplified  linearized  models  in  creeping 
flows  but  even  then  tedious  mathematical  manipulations  were  involved. 
Therefore  to  circumvent  this  problem,  we  assumed  that  the  length  of 
the  third  axis  was  given  by  the  geometric  mean  of  the  other  two 
measured  lengths.  In  brief,  the  volume  of  the  ellipsoid  was  taken  to 
be  proportional  to  the  cube  of  the  geometric  mean  radius,  1^.  Hence 
an  equivalent  expression  of  equation  (7)  for  the  ellipsoidal  bubbles 
is: 


P  <t)  -  <P0  *  P)  (VY  .  |S  .  (3) 

x  o  m  m 

Before  we  leave  this  topic,  we  would  like  to  point  out  that  the  error 
in  the  measurement  of  lengths  associated  with  the  bubbles  was  less  than 
2  urn. 

It  has  always  been  difficult  to  decide  on  the  value  of  the  poly¬ 
tropic  constant  y  to  be  used.  For  convenience,  many  researchers,  for 
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example,  Franklin  and  McMillan  (36,37) t  piesset  and  Prosperetti 
have  traditionally  taken  y  to  be  4/3.  This  value  was  also  used  in  the 
present  computation.  In  order  to  have  an  idea  on  the  lower  bound  of 
the  calculated  results,  the  isothermal  case  was  also  considered. 

Before  this  chapter  comes  to  a  close,  it  should  be  mentioned 
that  in  all  the  holograms  that  were  analyzed,  no  fragnentation  of  the 
bubbles  was  seen.  However,  a  few  badly  sheared  bubbles  were  observed, 
as  illustrated  in  Figure  4.6,  but  they  occur  very  infrequently.  The 
volumes  of  these  sheared  bubbles  were  difficult  to  determine  accurately 
and  were  omitted  in  the  computation.  It  has  also  been  discovered  that 
the  bigger  bubbles  have  a  greater  tendency  of  being  badly  sheared  than 
the  smaller  ones. 

Results  obtained  for  the  pressure  fluctuations  are  presented  and 
discussed  in  the  next  two  chapters. 
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5.  PRESENTATION  OF  RESULTS 


5.1  Introduction 

The  results  of  the  measured  cavitation  indices  and  flow  visuali¬ 
zation  are  presented  in  this  chapter,  Also  included  are  the  nuclei 
population  and  pressure  fluctuation  measurements  in  the  three  jets 
investigated.  The  main  features  of  the  various  plots  and  figures  are 
pointed  out.  However,  no  attempt  is  made  to  explain  these  features 
until  the  next  chapter  where  a  complete  discussion  is  presented. 

5.2  Measured  Scale  Effects 

The  results  of  the  measured  scale  effects  under  the  various  flow 
conditions  were  plotted  with  the  incipient  cavitation  index  as  the 
abscissa  and  the  jet  exit  Reynolds  number  as  the  ordinate.  These 
are  shown  in  Figures  5.1  to  5.24, 

Figures  5.1  to  5,4  show  the  effect  of  jet  size  at  different  dis¬ 
solved  air  content, a,  in  the  water  when  the  flow  is  unseeded.  An  im¬ 
mediate  observation  is  that  for  a  particular  nozzle  size,  the  inci¬ 
pient  index  is  independent  of  the  Reynolds  number.  It  can  also  be 
seen  that  while  the  two  smaller  jets  differ  only  slightly  in  their 
measured  indices,  at  the  same  degree  of  air  saturation,  the  largest  jet 
exhibits  an  incipient  index  which  is  significantly  higher  than  the 
other  two.  These  observations  are  true  over  the  entire  range  of  dis¬ 
solved  air  contents  used  presently.  It  should  be  pointed  out  that  the 
data  points  on  the  graphs  are  the  average  measured  value  and  the  bars 
denote  the  scatter.  Furthermore,  the  reason  why  relatively  few  data 
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points  were  obtained  for  the  cases  when  the  air  content  was  less  than 
8  ppm  was  because  the  jet  could  not  he  made  to  cavitate  at  the  lower 
velocities. 

The  above  data  are  replotted  as  shown  in  Figures  5.5  to  5,7  to 
demonstrate  the  effect  of  the  dissolved  air  content  on  the  susceptibility 
of  the  jet  to  cavitation  inception,  for  a  constant  nozzle  size.  These 
plots  clearly  illustrate  the  suppression  of  the  inception  index  with 
decreasing  air  content.  This  effect  is  much  more  pronounced  for  the 
3.17  mm  and  4.76  mm  jets  than  for  the  6.35  mm  jet. 

The  results  for  the  above  series  of  test  are  summarized  in 
Table  5.1  below  to  provide  a  convenient  reference. 


Size 

of 

Jet 

(Ml 

Measured  Value  of  Incipient  Cavitation 
Dissolved  Air  Content  (/ppm) 

14.1  :  10. 9  7.7 

Index 

4.2 

Comments 

3.17 

0.16 

0.08 

0.06 

0.03 

Independent 

4.76 

0.15 

0.11 

0.07 

0,06 

of 

Jet  Exit 

6.35 

Q.27 

0.24 

0.23 

0.22 

Reynolds 

Number 

Table  5.1.  Summary  of  the  measured  cavitation  indices  for 
cavitating  jets. 


One  objective  of  the  present  study  was  to  determine  the  effect 
of  nuclei  population  on  cavitation  inception.  This  was  accomplished 
by  seeding  the  flow  with  electrolysis  bubbles.  Electrolysis  currents 
ranging  from  0,05  to  5  mA  were  used  to  artificially  create  different 
nuclei  populations.  The  results  of  this  study  are  presented  in  Figures 
5.8  to  5,10.  For  comparison,  the  inception  data  obtained  at  the  same 


dissolved  air  content  but  with  no  flow  seeding  are  also  included  in 
the  plots,  rt  can  Be  seen  from  these  plots  that  the  inception  number, 
<y.j ,  for  both  the  4.76  mm  (figure  5,9i  and  6,35  mm  jets  (figure  5,10) 
exhibit  a  similar  behavior.  This  is  to  say  that  for  a  constant  nozzle 
size,  the  a.  display  an  increasing  dependence  on  Reynolds  number 
with  increasing  electrolysis  current,  r.  This  dependence  of  cr^  on 
Reynolds  number  does  not  go  on  indefinitely  with  increasing  I,  but 
approaches  a  limit  when  the  current  exceeds  a  certain  value.  This 
limit  is  reached  when  I  >  0,4  mA  for  the  4.76  mm  jet,  and  1  >  0,6  mA 
for  the  largest  jet.  As  for  the  smallest  jet,  it  can  be  seen  from 
Figure  5.8(a)  that  the  behavior  described  above  for  the  other  two  jets 
is  duplicated  here  at  an  air  content  of  14  ppm,  with  the  Reynolds 
number  dependence  flattening  out  at  I  >  0.2  mA.  At  the  other  three 
air  contents  (Figures  5.8(b)  to  5.8(d)),  the  Inception  number  of  the 
3.17  mm  jet  does  not  show  any  significant  dependence  on  Reynolds 
number,  over  the  range  of  electrolysis  currents  used.  However,  the 
measured  Incipient  index  increases  with  increasing  current.  At  this 
juncture.  It  should  be  pointed  out  that  the  scatter  on  the  measured 
inception  number  for  all  the  tests  with  flow  seeding  is  approximately 
±  0.02.  It  is  also  worth  noting  that  the  experimental  data  for  a 
particular  nozzle  size  were  measured  over  as  wide  a  range  of  jet 
velocities  (or  Reynolds  nunber)  as  possible;  the  limit  was  imposed  by 
the  operating  range  of  pressures  in  the  LTWT  which  was  approximately 
0.2  to  1.0  atmospheres. 
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The  experimental  data  are  replotted  as  shown  in  Figures  5,11  to 
5,16  to  elucidate  the  effect  of  the  dissolved  air  content  on  at  a 
constant  "seeding  rate".  These  plots  clearly  demonstrate  that  when 
the  flow  Is  seeded,  the  influence  on  cr^  is  determined  by  the  electrolysis 
current  I,  and  does  not  depend  on  the  dissolved  air  content.  Excep¬ 
tions  to  this  observation  occur  for  the  largest  jet  at  I  £  0.6  mA 
CFigures  5.15(a)  to  5 . 1 5 Cd i  1  and  at  I  =  5.0  mA  (Figure  5.16(d))..  In 
these  instances,  the  amount  of  dissolved  air  in  the  jet  water  appear 
to  have  some  effect  on  the  measured  inception  number,  although  no 
regular  pattern  could  he  observed. 

Figures  5.17  to  5,24  were  plotted  with  the  Intention  of  presenting 
the  experimental  results  in  such  a  way  that  the  influence  of  jet  size 
on  cr.  for  jets  operating  under  similar  conditions  (some  dissolved  air 
content  and  electrolysis  current)  could  be  easily  grasped.  On 
examination  of  Figures  5.17  and  5.18,  it  can  be  seen  that  at  an  air 
content  of  approximately  14  ppm  and  low  values  of  I,  the  inception 
data  for  the  smallest  jet  (3.17  mm)  appear  to  serve  as  an  extension 
to  the  results  of  the  4,76  mm  jet  In  the  lower  Reynolds  number  range. 
However,  as  I  is  Increased,  these  data  begin  to  deviate  from  the  data 
of  the  4.76  mm  jet  and  begin  to  approach  those  of  the  largest  jet  at 
lower  Reynolds  number. 

At  an  air  content  of  IQ, 3  ppm,  the  inception  data  for  the  two 
smaller  jets  again  show  signs  of  supplementing  each  other  over  the 
range  of  ReD  Investigated,  as  shown  graphically  in  Figures  5.19  and 
5.20.  This  observation  is  especially  true  for  I  <  0.6  mA.  At  the 
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same  time,  the  largest  jet  exhibits  a  hiS^r  y$lue  of  than  the 
two  smaller  jets  at  corresponding  Reynolds  number,  It  is  also 
Interesting  to  note  that  the  slope  of  the  experimental  results  for 
the  two  larger  nozzles  are  nearly  the  same  in  the  range  of  electrolysis 
current  between  0,05  mA  to  2,0  mA  (figures  5,19(al  to  5,2Q(c).l. 

On  close  inspection  of  Figures  5,21  to  5,24,  it  can  be  in¬ 
ferred  that  at  both  air  content  of  7,7  ppm  and  4,2  ppm,  there  are 
instances  when  the  data  of  the  two  smaller  jets  supplement  each 
other.  At  other  instances,  the  jets  portray  different  behavior.  In 
any  event,  the  measured  for  the  largest  nozzle  is  always  higher 
than  for  the  other  two,  at  the  same  Reynolds  number. 

Tn  addition  to  the  observations  described  above,  flow  visualiza¬ 
tion,  nuclei  population  density  and  pressure  fluctuation  measurements 
in  the  jets  were  also  carried  out.  These  were  made  in  the  hope  that 
the  information  obtained  would  provide  some  probable  explanations 
for  the  observed  scale  effects.  The  results  for  these  tests  are 
presented  in  the  next  three  sections  of  this  chapter, 

5.3  Flow  Visualization 

5.3,1  Observation  under  Stroboscopic  Lighting 

As  the  jet  was  made  to  cavitate,  the  following  sequence  of 
events  was  observed  under  stroboscopic  lighting. 

Small  irregular  shaped  cavities  were  seen  in  the  jet 
at  inception.  Schlieren  photographs  corresponding  to  this 
Situation  are  shown  in  Figures  5.25,  5.26  and  5.27 
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where  the  cayities  appear  in  the  form  of  black,  unevenly  shaped 
spots.  It  should  Be  pointed  out  that  the  direction  of  flow  of  the 
jet  in  the  above-mentioned  figures  Is  from  left  to  right.  These 
cavities  were  transient  In  nature  and  they  collapsed  and  new  ones 
were  formed  randomly  in  the  jet,  within  an  axial  range  of  6  to  16 
jet  diameters  downstream  of  the  nozzle  lip.  Tn  some  runs,  cavities 
were  seen  to  appear  as  near  as  four  jet  diameters  and  as  far  as 
twenty  diameters  downstream  but  their  occurrences  were  relatively 
infrequent.  Inception  was  accompanied  by  "intermittent  crepitant 
sound"  caused  by  the  random  exploding  and  collapsing  cavities.  The 
"Intermittent  crepitant  sound"  continued  for  a  short  duration,  but 
with  Increasing  frequency,  as  the  pressure  in  the  tunnel  was  reduced 
beyond  the  Inception  point,  until  an  almost  "continuous  popping 
sound"  was  heard.  Such  a  situation  is  recorded  in  Figure  5.28.  The 
cavities  were  now  more  extensive  and  much  larger  in  size  than  at 
inception. 

With  further  reduction  in  pressure,  the  almost  "continuous 
popping  sound"  gave  way  to  a  steady  roar.  The  jet  was  then  said  to 
be  fully  cavltatlng  since  no  more  changes  in  the  appearance  of  the 
cavltatlng  jet  were  observed  beyond  this  point.  At  this  stage, 
almost  the  entire  region  of  the  jet  between  one  to  twenty  diameters 
downstream  were  covered  by  cavities.  These  cavities  have  been 
observed  to  collapse  Into  a  cloud  of  very  fine  "mist"  at  distances 
greater  than  about  twenty  diameters.  All  these  features  are  clearly 
depicted  in  Figures  5.29(a)  and  0>). 
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The  sequence  of  events  described  above  was  true  when  the  jet  was 
nearly  saturated  with  air  or  when  the  flow  was  seeded.  However,  at  an 
air  content  of  less  than  9  ppm  and  with  no  flow  seeding,  the  cavita¬ 
tion  process  was  quite  different.  Inception  occurred  in  a  violent 
burst  of  cavities,  as  recorded  in  Figures  5.30(a),  (b)  and  (c). 

A  comparison  with  the  corresponding  cases  of  the  high  air  content 
jets  revealed  that  these  cavities  were  more  extensive.  It  was 
observed  that  the  cavities  normally  lasted  for  a  very  short  time,  and 
collapsed  soon  after  formation.  If  the  pressure  in  the  water  tunnel 
was  kept  at  the  Inception  condition  for  a  prolonged  period  or  lowered 
past  this  point,  a  second  and  subsequent  bursts  of  cavities  would 
result.  There  appeared  to  be  no  regular  time  interval  between  bursts, 
when  the  jet  was  effectively  a  single  phase  water  jet. 

On  further  reduction  in  pressure,  an  instant  arose  when  the  jet 
became  fully  cavltatlng.  This  was  seen  only  when  Vj  >  25  m/s.  (.for  the 
two  smaller  jets).  However,  there  was  a  difference  between  this  case 
and  the  high  air  content  jet  described  earlier.  The  low  air  content 
jet  did  not  remain  fully  cavitating  over  a  continuous  period  but 
alternated  between  a  fully  cavitating  and  a  non-cavitating  single 
phase  jet.  Again  these  alternating  actions  of  the  jet  occurred  at 
Irregular  Intervals. 

All  the  above  description  was  true  of  the  sequence  of  events 
that  was  exhibited  by  the  three  jets  under  Investigation,  with  the 
exception  that  In  the  6.35  mm  jet,  cavities  at  inception  were  concen¬ 
trated  in  the  region  of  two  to  fourteen  diameters  from  the  nozzle  exit. 
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It  should  ail  so  he  mentioned  that  In  some  test  runs  at  a  dis¬ 
solved  air  content  of  4  ppm,  the  two  smaller  jets  could  not  be  made 
to  cavitate  at  all,  Below  a  jet  exit  velocity  of  27  m/s,  even  when 
the  jets  were  kept  at  the  lowest  attainable  tunnel  pressure  over  a 
period  of  approximately  sixty  seconds.  However,  when  the  flow  was 
seeded  by  passing  a  very  small  amount  of  current  C0.Q1  mA}  through 
the  electrode,  cavitation  in  the  jet  occurred.  The  cavitation  ceased 
almost  instantly  when  the  current  was  turned  off.  This  clearly 
demonstrated  the  role  of  nuclei  in  the  inception  process. 

5,3.2,  Flow  Observation  with  Schlieren  Photography 

Schlieren  photography  was  utilized  as  a  means  to  give  a 
permanent  record  of  the  appearance  of  the  cavitating  and  non-cavitating 
jets.  The  photographic  records  of  the  jets  at  and  beyond  the  inception 
point  confirmed  the  observations  made  under  stroboscopic  lighting. 

As  mentioned  in  Section  3,4,  schlieren  photographs  were  taken 
with  the  knife  edge  at  two  different  positions.  This  was  done  to 
ensure  that  all  the  features  in  the  jet  were  revealed.  Figures 
5,31  Cal  -  (d)  are  a  series  of  schlieren  photographs  of  the  *.35  mm 
CO. 25  inch)  jet  at  different  exit  velocities,  taken  with  the  knife 
edge  parallel  to  the  direction  of  flow.  An  obvious  observation  Is 
that  the  jets  are  laminar  at  exit  and  the  length  of  the  laminar  region 
retreats  towards  the  nozzle  lip  as  the  velocity  in  the  jet  Increases. 
Another  striking  feature  of  the  jets  is  the  appearance  of  cellular 
shaped  structures  at  the  end  of  the  laminar  region  just  upstream  of 
the  turbulent  regime.  On  closer  Inspection,  it  can  be  seen  that  the  size 
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of  these  structures  decreases  with  increasing  velocity, 

Figures  5,32Cal  -  (di  are  of  the  same  jet  But  recorded  with  the 
cut-off  edge  perpendicular  to  the  flow  direction.  All  the  above  des¬ 
cribed  features  are  again  seen  here  except  that  in  the  turbulent 
region  of  the  jet,  the  scale  of  the  grainy  structures  is  a  little 
finer.  In  addition,  the  jet  taken  with  this  knife  edge  position  has 
a  "flat"  appearance,  whereas  in  the  earlier  situation  it  appears  more 
"tubular"  in  form. 

Corresponding  sets  of  figures  are  also  presented  for  the  4.76 
nin  CO -188  inch)  and  3,17  mm  CO. 125  inch)  jets.  These  are  found  in 
Figures  5.33Ca)  -  Cd}  and  5.34Cal  -  Cd)  for  the  4.76  mm  jet  and 
Figures  5,35Ca)  -  Cd)  and  5 , 36 Ca 1  -  Cd)  for  the  smallest  jet.  A 
careful  examination  of  these  figures  shows  that  the  features  exhibited 
by  the  6,35  mm  jet  are  also  displayed  by  these  two  smaller  jets,  but 
are  correspondingly  smaller  in  scale. 

Two  types  of  measurements  were  obtained  from  the  schlieren  photo¬ 
graphs.  The  first  of  these  measurements  was  the  determination  of  the 
non-dimensional  length  of  the  laminar  portion, L./D,  of  the  jet  as  a 

J 

function  of  the  jet  exit  Reynolds  number.  These  results  are  summarized 
graphically  in  Figures  5,37Ca),  (b)  and  Cc).  It  can  be  inferred  from 
these  plots  that  for  the  two  larger  jets,  Lj/D  decays  exponentially 
with  increasing  Reynolds  number  while  the  smallest  jet  shows  a  linear 
dependence,  at  least  over  the  range  of  Red  currently  investigated. 
Secondly,  the  spread  rate  of  the  jets  was  determined.  This  was  done, 
for  each  jet,  by  measuring  the  width  of  the  jet  at  various  down- 
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stream  stations  in  the  jet,  from  the  photographic  records,  The  results 
were  plotted  in  a  non-dimensional  form  as  shown  in  Figures  5.38(a).  - 
(cl.  It  is  evident  from  these  graphs  that  the  spread  of  the  two 
smaller  jets  is  Independent  of  the  jet  velocity.  However,  there  are 
indications  that  the  spread  of  the  largest  jet,  especially  near  the 
nozzle  exit,  is  affected  by  the  exit  velocity.  It  should  be  mentioned 
that  each  data  point  on  the  plots  corresponded  to  the  average  of 
three  measurements  made  from  identical  runs. 

5,3.3  Flow  Visualization  with  Holography 

Holographic  studies  of  the  non-cavitating  jet  at  various 
cavitation  number  revealed  that,  irrespective  of  the  dissolved  air 
content  in  the  water,  no  nuclei  greater  than  40  ym  in  diameter  were 
found  in  any  region  of  the  three  jets  investigated.  This  observation 
was  even  true  at  flow  conditions  approaching  the  cavitation  inception 
point.  It  should  be  noted  that  the  size  of  the  holograms  permitted 
only  5  cm  of  the  downstream  distance  of  the  jet  to  be  scrutinized. 

When  the  flow  was  seeded,  a  large  amount  of  bubbles  of  approxi¬ 
mately  10  urn  in  diameter  were  detected  in  the  jet.  Apparently,  the 
size  of  these  bubbles  was  not  significantly  affected  by  either  the  jet 
velocity  or  the  magnitude  of  the  electrolysis  current. 

A  careful  examination  of  holograms  of  jets  at  inception  con¬ 
firmed  that  the  cavities  were  contained  In  the  same  distance  range 
mentioned  In  earlier  sections  of  this  chapter.  However,  it  was  also 
found  that  the  cavities  could  occur  at  any  region  across  the  width 
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of  the  jet  favoring  no  particular  radial  position,  The  aboye  findings 
were  obtained  through  the  analysis  of  18  holograms. 


5.4.  Nuclei  Population  Density  Measurements  by  Holography 

The  nuclei  populations  in  the  jets  were  measured  from  holo¬ 
grams  of  the  non-cavitating  jet  at  different  values  of  a  and  at 
four  different  dissolved  air  content.  In  a  typical  hologram  where  the 
nuclei  population  was  to  be  determined,  the  jet  was  divided  into  four 
different  regions  labelled  1,  2,  3  and  4  as  shown  schematically  in 
Figure  5.39,  The  nuclei  count  in  each  of  these  regions  was  reduced 
to  a  number  density  distribution  function  N(R1  defined  as  follows; 


NOD* 


number  of  nuclei  per  unit  volume  with  radii  between  and  R2 
- - 


where  the  mean  bubble  radius,  R,  is  given  by  (R-j  +  R2)/2.  Note  that  R1 
and  R2  correspond  to  the  lower  and  upper  limits  of  a  certain  size 
range.  The  nuclei  were  sorted  into  the  following  categories: 

5  pm  <  R  10  ym,  10  ym  <  R  <_  20  pm  and  R  >  20  ym.  Although  nuclei 
less  than  5  ym  in  radius  were  observed,  they  were  not  recorded  in  the 
present  work,  because  they  could  not  be  counted  with  the  same  degree 
of  certainty  as  the  bigger  size  nuclei  due  to  optical  noise.  In  all 
the  nuclei  measurements  made  presently,  the  control  volume  that  was 
used  was  defined  as 


Control  Volume  sDxDxO 
where  0  is  the  jet  size. 
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Figures  5,4Q(_al  -  (dl,  5,41  Cal  -  (dl  and  5,42(al  -  CpL  are  plots 
of  the  num6er  density  distribution  of  the  3,17  mm,  4,76  mm  and  the 
6.35  mm  jet,  respectively.  The  plots  in  each  set  represent  the 
density  functions  in  the  designated  regions  of  the  jet  and  are 
labelled  accordingly  on  the  graphs.  It  can  Be  seen  from  these  plots 
that  there  is  a  definite  dependence  of  the  numBer  density  function 
on  the  dissolved  air  content, a.  Also,  it  can  Be  inferred  from  the 
figures  that  the  data  points  corresponding  to  a  specific  air  content 
do  not  show  a  regular  pattern  with  which  the  measured  N(Rfl  varies 
with  the  cavitation  number.  In  addition,  a  comparison  of  the  data 
from  plots  of  similarly  labelled  regions  in  the  three  jets  indicated 
that  at  a  constant  air  content,  the  measured  values  of  NCR)  for  a 
particular  size  category  fall  approximately  in  the  same  range  of 
values.  Consequently,  the  value  of  NCR).  for  a  specific  air  content 
was  obtained  by  averaging  the  results  of  all  the  jets  at  that  air 
content.  In  so  doing,  it  was  assumed  that  the  density  function  did 
not  depend  on  the  cavitation  number.  This  appears  to  be  a  fair 
assumption  since  as  pointed  out  earlier,  the  dependence  of  N(J)  on  a 
does  not  follow  a  uniform  pattern.  It  should  be  mentioned  that 
documentation  of  the  nuclei  population  in  Region  4  of  the  largest  jet 
was  not  possible  because  that  region  was  not  contained  in  the  holo¬ 
graphic  record,  due  to  the  small  size  of  the  hologram. 

Table  5.2  summarizes  the  results  of  the  present  nuclei  popula¬ 
tion  density  distribution  measurements.  The  entries  in  this  tabulation 
will  be  analyzed  and  discussed  in  the  next  chapter.  For  now,  it  suf- 
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Table  5.2.  Summary  of  the  number  density  function  of  nuclei  measurements  In  jets. 
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flees  to  say  that  no  nucleus  with,  a  radius  greater  than  2Q  pm  was 
sighted  in  any  of  the  jets. 

When  electrolysis  6u6bles  were  introduced  into  the  flow,  it 
was  found  that  the  nuclei  were  so  numerous  that  the  number  could  not 
be  counted  with  any  degree  of  certainty.  Hence  the  effort  was 
abandoned.  However,  an  attempt  was  made  to  estimate  the  number  of 
bubbles  that  were  generated  by  the  electrolysis  current.  This  estimate 
was  based  on  the  work  of  Kufper  and  is  considered  in  Chapter  6. 

5.5  Measurements  of  Pressure  Fluctuation 

The  motivation  behind  the  desire  to  measure  the  pressure  fluc¬ 
tuations  in  a  jet  and  the  techniques  employed  to  achieve  this  objective 
have  already  been  discussed  in  Chapter  4,  The  results  are  presented 
below. 

Figures  5,43  and  5.44  are  examples  of  the  instantaneous  pres¬ 
sure  fluctuations  in  the  3.17  irm  and  4,76  mm  jets,  respectively.  The 
data  for  the  6.35  mm  jet  are  not  presented  because,  as  pointed  out 
earlier,  the  • ize  of  the  hologram  prevented  any  detail  mapping  of 
the  pressure  field.  In  these  graphs,  the  pressure  fluctuation  P' 

Is  expressed  as  a  fraction  of  the  dynamic  head  at  the  exit  of  the 
jet  and  is  plotted  against  the  non-dimensiional  axial  position.  The 
lines  on  the  graphs  join  data  points  at  constant  radial  position  in 
the  jet.  The  plots  were  made  on  a  predetermined  radial  grid  spacing 
as  shown  on  the  graph  and  all  the  data  points  that  fell  within  a 
value  of  ±0,1  of  a  particular  non-dimensional  radial  position  were 
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assumed  to  he  at  that  position,  These  graphs  clearly  show  that  a 
whole  spectrum  of  pressure  fluctuations  ranging  from  0  percent  to  16Q 
percent  of  the  dynamic  head  Is  present  at  any  one  tnstant  in  the  jet. 
Figure  5,45  is  a  proBability  histogram  which  shows  the  distri¬ 
bution  of  the  Intensities  of  pressure  fluctuations  in  the  jet.  This 
curve  was  obtained  by  incorporating  all  the  data  from  17  holograms  and 

was  plotted  in  a  non-dimensional  form  as  shown,  being  the  number 

★ 

of  data  points  in  a  certain  category  and  (=130)  corresponds  to  the 
maximum  number  recorded  on  the  probability  histogram.  It  should  be 
pointed  out  that  only  the  data  on  pressure  fluctuations  outside  of 
the  potential  core  region  of  the  jet  were  used  in  the  present  plot. 
This  was  because  within  the  potential  core  little  or  no  pressure 
fluctuations  were  registered.  Hence,  if  these  data  points  were  in¬ 
cluded  in  the  histogram,  they  would  have  caused  an  unfair  bias  towards 

the  lowest  pressure  fluctuation  intensity  category.  It  can  be  in¬ 
ferred  from  this  graph  that  the  pressure  fluctuation  has  an  almost 
Gaussian  distribution,  but  is  skewed  somewhat.  Notice  also  that  the 
fluctuations  peak  at  a  higher  positive  value  than  a  negative  one. 
Figure  5.46  is  the  distribution  curve  when  the  data  are  replotted  as 
a  function  of  the  modulus  of  the  dimensionless  pressure  fluctuation 
intensity.  Both  Figures  5,45  and  5,46  are  for  the  case  when  the 
pressure  calculations  were  made  with  a  polytropic  index  of  4/3,  For 
comparison, corresponding  plots  for  the  isothermal  case  are  found  in 
Figures  5.47  and  5,48, 
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Fig.  5.2.  Variation  of  the  Incipient  cavitation  Index  with  jet  exit  Reynolds  number  for  the  3.17mm, 
4.76mm  and  6.35mm  jet  at  a  dissolved  air  content  of  10.9ppm. 
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Fig.  5.3.  Variation  of  the  Incipient  cavitation  Index  with  jet  exit  Reynolds  number  for  the  3.17mm 
4.76mm  and  6.35mm  jet  at  a  dissolved  air  content  of  7.7 ppm. 
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Fig.  5.4.  Variation  of  the  incipient  cavitation  index  with  jet  exit  Reynolds  number  for  t 
and  6.35mm  jet  at  a  dissolved  air  content  of  4.3ppm. 
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Fig.  5.5.  Summary  of  the  variation  of  the  incipient  cavitation  index  with  Reynolds  number  at  different 
dissolved  air  contents  for  the  3.17mm  jet. 
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Fig.  5.6.  Summary  of  the  variation  of  the  incipient  cavitation  index  with  Reynolds  number  at  different 
dissolved  air  contents  for  the  4.76mm  jet. 
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Fit.  5.7.  Summary  of  the  variation  of  the  incipient  cavitation  index  with  Reynolds  number  at  different 
dissolved  air  contents  for  the  6.35mm  jet. 


Fig.  5.8.  Variation  of  the  incipient  cavitation  index  with  Reynolds  number  at  different  electrolysis  currents 
for  the  3.17mm  jet:  (a)  a  =  1 3 . 98ppir. ,  (b)  a=11.01ppm,  (c)  a  =  7.53ppm,  and  (d)  a  =  4.23ppm. 


Variation  of  the  Incipient  cavitation  Index  with  Reynolds  number  at  different  electrolsls  currents 
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Illustrating  the  effect  of  dissolved  air  content  on  the  Incipient 
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Fig.  5.19.  Graph  of  oj  versus  Ren  showing  the  Influence  of  Jet  size  on  the  Incipient  cavitation  Index  at 
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Fig.  5.28.  Schlieren  photograph  taken  at  the  onset  of  the  almost 
"continuous  popping  sound"  state  of  cavity  development 
the  6.35mm  jet  at  Re^  1 . 29  x 105  and  o =  0.19. 


Fig.  5.29(b).  Schlieren  photograph  of  the  same  jet  as  in  (a)  above 
showing  the  region  where  the  cavity  collapsed. 


Fig.  5.30.  Schlieren  photoqraphs  of  cavitation  inception 
jet  at  Ren  =  2.06x10  and  =  0.22,  (b)  4. 76m" 
ReD  =  1.47  jclO'  and  ^  =  0.09,  (c)  3.17mm  jot  a 
and  =0.08.  The  dissolved  air  content  is  7. 
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Fig.  5.37.  Graphs  of  U/D  vs  Rep  for  the  (a)  6.35mm  jet,  (b)  4.76mm  jet 
and  (c)  3.17mm  jet. 
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Fig.  5.38.  Plots  of  the  non-dimensional  jet  width  6/D  versus  the  di¬ 
mensionless  axial  position  for  the  (a)  3.17mm  jet,  (b)  4.76irm 
jet,  and  (c)  6.35mm  jet. 
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Fig.  5.39.  A  schematic  drawing  depicting  the  4  regions  in  the  jet  where  nuclei  counts  were  carried  out 
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Fig.  5.40.  Plots  of  the  nuclei  distributions  for  the  3.17mm  jet: 

(a)  Region  1,  (b)  Region  2,  (c)  Region  3,  and  (d)  Region  4 
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Fig.  5.44.  Graphs  of  the  instantaneous  local  pressure  fluctuations  as  a  function  of  the  dimensionless 
axial  position  in  the  4.76mm  jet. 
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Probability  histogram  showing  the  distribution  of  the  pressure  fluctuation  intensities 
(Polytropic  index  of  4/3) 
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Fig.  5.46.  Probability  histogram  showing  the  distribution  of  the  modulus  of  pressure  fluctuation  intensities 
in  a  jet.  (Polytropic  index  of  4/3). 
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Probability  histogram  showing  the  distribution  of  the  pressure  fluctuation  intensities 
in  a  jet  (Isothermal  case). 
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Fig.  5.48.  Probability  histogram  showing  the  distribution  of  the  modulus  of  pressure  fluctuation  intensities 
in  a  jet.  (Isothermal  case). 
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6.  DISCUSSION 


6.1  Introduction 

The  results  of  the  present  studies  are  quantitatively  and  quali¬ 
tatively  discussed  here.  Explanations  for  the  causes  of  certain  ob¬ 
served  scale  effects  are  provided  through  correlation  of  all  the  test 
results.  In  situations  where  insufficient  data  are  available  to  make 
definite  statements,  probable  reasons  or  explanations  are  tendered 
instead. 

This  chapter  is  divided  into  four  sections.  The  first  section 
deals  with  the  results  of  the  unseeded  jets  while  the  second  focuses 
on  the  jets  when  flow  seeding  was  carried  out.  This  is  followed  by  a 
comparison  of  results  with  the  work  of  others.  Finally,  the  chapter 
closes  with  proposals  for  future  work. 

6.2  Unseeded  Jets 

Results  presented  in  Table  5.1  of  Chapter  5  have  shown  that  for 
a  jet  of  constant  size,  the  incipient  cavitation  index  decreases  with 
decreasing  dissolved  air  content  in  the  water.  It  has  long  been  postu¬ 
lated  that  this  suppression  of  the  inception  point  is  due  to  a  reduc¬ 
tion  in  the  number  of  cavitation  nuclei.  However,  this  has  never  been 
proven  or  even  documented  quantitatively  prior  to  the  present  work, 
insofar  as  studies  on  jets  are  concerned.  The  tabulation  of  the  nuclei 
population  density  (Table  5,2}  clearly  shows  that  when  the  dissolved 
air  content  is  reduced  from  14.3  ppm  to  4.2  ppm,  the  nuclei  population 
is  reduced  by  a  factor  of  approximately  35,  thus  supporting  the  postu- 
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late,  On  closer  inspection  of  the  entries  In  Table  5,2,  it  can  also 
be  seen  that  the  nuclei  population  density  is  higher  in  regions  out¬ 
side  the  jet  CReglons  2  and  41  than  those  inside  it.  This  may  seem  a 
little  puzzling  at  first,  but  as  mentioned  earlier  the  jet  in  the  pre¬ 
sent  studies  was  created  by  discharging  a  column  of  water  from  a  pres¬ 
surized  vessel.  Hence  it  is  to  be  expected  that  some  of  the  gas 
bubbles  in  the  jet  were  forced  into  solution  by  the  pressurization 
process.  Consequently  the  nuclei  population  near  the  jet  exit  (Region 
1)  would  be  lower  than  in  the  entrainment  zone  (Region  2).  One 
would  also  expect  the  nuclei  count  to  increase  with  downstream  dis¬ 
tance  from  the  nozzle  lip  because  of  flow  entrainment  from  the  sur¬ 
rounding  fluid  which  contains  a  higher  nuclei  population.  This  can 
be  shown  to  be  the  case  by  comparing  the  nuclei  population  density  in 
Region  3  to  Region  1.  It  is  worth  pointing  out  that  as  the  dissolved 
air  content  drops  so  does  the  difference  in  the  nuclei  population  in 
regions  inside  and  outside  the  jet. 

It  has  also  been  found  in  the  same  nuclei  population  measurements 
discussed  above  that  the  population  density  of  a  fixed  dissolved  air 
content  is  independent  of  the  cavitation  index  in  the  jet,  at  least 
up  to  the  inception  point.  It  is  believed  that  this  finding  is  not 
universally  true  but  is  a  consequence  of  the  test  procedures  cur¬ 
rently  practiced.  Remember  that  the  jet  was  discharged  into  a 
stationary  tunnel  sitting  at  an  absolute  pressure  of  approximately  13 
psia  and  made  to  cavitate  by  rapidly  lowering  the  pressure  in  the  tunnel. 
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Since  the  duration  of  a  typical  run  was  less  than  10  seconds,  it  is 
fair  to  assume  the  pressure  reduction  process  itself  would  not  sig¬ 
nificantly  affect  the  nuclei  population,  rt  is  our  opinion  that  the 
nuclei  count  at  a  particular  air  content  under  such  an  operating 
condition  is  determined  by  the  steady  state  condition  in  the  water 
tunnel,  namely,  the  pressure  at  which  the  water  tunnel  was  kept  prior 
to  a  test  run. 

The  results  as  tabulated  in  Table  5.1  in  Section  5.2  show  that 
the  measured  incipient  cavitation  indices  for  the  3.17  mm  and  4.76  mm 
jets  at  the  same  dissolved  air  content  are  not  too  different  from 
each  other.  On  the  other  hand,  the  6.35  mm  jet  exhibits  a  cavitation 
inception  number  that  is  substantially  different  from  the  other  two. 

The  most  probable  explanation  for  this  observation  would  be  the  dif¬ 
ferent  nuclei  population  in  the  cavitation-prone  region  of  the  jet. 
However,  to  account  for  the  big  difference  in  the  measured  a^'s,  the 
nuclei  population  in  the  largest  jet  should  be  significantly  larger 
than  the  other  two.  An  attempt  will  now  be  undertaken  to  demonstrate 
that  this  is  the  case. 

It  has  already  been  established  that  at  any  Jissolved  air  content, 
the  nuclei  population  density  is  independent  of  the  size  of  the  jet. 

It  follows  then  that  at  any  instant  the  number  of  nuclei  present  in 
the  region  that  is  susceptible  to  cavitation  inception  is  proportional 
to  the  volume  of  liquid  in  that  region.  As  has  been  reported  earlier, 
the  cavitation-prone  region  of  the  two  smaller  jets  extends  from  6  to 
16  jet  diameters  downstream  of  the  nozzle  exit  while  in  the  largest  jet 
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it  covers  a  distance  of  2  to  14  diameters,  Using  these  results  and 

the  data  on  the  spread  of  the  three  jets  as  presented  in  Figures  5.38(a)-(c), 

it  was  calculated  that  the  relevant  volumes  in  the  6.35  mm 

and  4,76  mm  jets  were  respectively  9.3  and  2,3  times  more  than  in  the 

smallest  jet.  The  volume  in  the  smallest  jet  was  computed  to  be  5.41 

3 

cm  .  Therefore,  for  the  reason  stated  at  the  outset  of  this  paragraph, 
it  follows  that  the  measured  cr^ ' s  are  directly  proportional  to  the 
number  of  nuclei  present. 

If  the  nuclei  population  was  the  only  determining  parameter,  then 
one  would  expect  the  measured  o\'s  to  be  scaled  up  or  down  in  the  same 
proportion  when  the  nuclei  population  is  altered  as  a  result  of  using 
water  with  different  dissolved  air  content.  However,  it  is  evident 
from  the  present  results  that  this  is  not  the  case;  the  effect  of  a 
change  in  the  number  of  nuclei  on  the  cavitation  inception  number  is 
more  pronounced  on  the  3.17  mm  and  4.76  mm  jets  than  for  the  6.35  mm 
jet.  Obviously  this  implies  that  at  least  one  other  scaling  para¬ 
meter  is  involved.  The  most  likely  candidate  is  the  pressure 
fluctuation  since  it  is  responsible  for  initiating  cavitation. 

Typically,  a  whole  spectrum  of  pressure  fluctuation  intensities 
exists  in  a  jet  CArndt  and  George  In  the  present  investigation, 

negative  peak  pressure  fluctuations  from  0  to  120  percent  of  the 
dynamic  head  of  the  jet  at  exit  were  detected.  However,  this  is  not 
to  say  that  cavitation  inception  would  occur  at  the  first  negative 
pressure  fluctuation  a  bubble  sees.  In  reality,  two  conditions  have 
to  be  present  before  inception  can  occur.  The  first  condition  is  that 
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the  pressure  fluctuation  must  be  intense  enough  to  create  a  local 
pressure  that  is  conducive  to  the  explosive  growth  of  the  bubble.  It 
can  be  shown  from  bubble  dynamics  that  this  local  pressure  is  reached 
when 

"  PV  ‘  W  (6J) 

The  second  condition  to  be  satisfied  is  that  the  local  pressure  P 
given  by  Equation  (6.1}  above  must  persist  at  or  even  below  that  value 
for  a  period  of  time  that  is  long  in  comparison  to  the  response  time 
of  the  bubble.  It  has  already  been  reported  that  no  nuclei  greater 
than  20  urn  in  radius  was  ever  detected  in  the  jets  currently  investi¬ 
gated.  Since  a  bubble  of  the  size  reported  has  a  response  time  of 
about  10  microseconds, which  is  about  two  orders  of  magnitude  shorter 
than  the  time  scale  of  the  pressure  fluctuations  or  the  traverse  time 
across  the  cavitation-prone  region,  we  would  expect  inception  to  take 
place  once  the  local  pressure  satisfies  Equation  (6.1).  It  is  also 
evident  that  the  pressure  fluctuations  recorded  in  the  present  tests 
are  intense  enough  to  account  for  cavitation  inception  in  all  the  three 
jets  studied. 

It  has  been  found  that  the  mean  static  pressure  at  which  cavitation 
inception  occurred  was  a  function  of  numerous  variables,  for  example, 
the  velocity  and  size  of  the  jet.  This  meant  that  different  minimum 
magnitudes  of  negative  peak  pressure  fluctuations  were  required  to 
initiate  cavitation  in  each  case.  However,  the  probability  histogram 
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on  the  pressure  fluctuations  in  Figure  5, 45  (a)  shows  that  the  different 
pressure  fluctuation  intensities  have  dissimilar  probabilities  of 
occurring.  Consequently,  the  cumulative  probability  of  a  favorable 
pressure  fluctuation  is  also  different  from  one  flow  condition  to 
another.  It  is  possible  that  this  variable  enters  as  a  second  scaling 
parameter. 

We  define  a  new  parameter  Np  as 

Np  =  n  x  (6,2) 

where  n  *  total  number  of  nuclei  present  in  the  region  of  the  jet  that 
is  susceptible  to  cavitation;  and  P(P^)  =  cumulative  probability  of 
favorable  negative  peak  pressure  fluctuation  intensities  that  would 
cause  the  explosive  growth  of  a  cavitation  nucleus. 

In  general,  nuclei  of  various  sizes  are  present  in  the  jets. 

These  nuclei  are  often  sorted  out  in  different  size  categories  to 
facilitate  documentation  of  their  presence.  In  addition  to  having 
different  population  densities,  the  nuclei  in  the  different  categories 
require  different  minimum  pressure  fluctuations  to  cause  them  to  grow 
explosively.  Consequently  the  value  of  Np  as  defined  in  Equation 
(6.2)  would  vary  from  one  size  category  to  another.  Hence,  Equation 
(6.2)  should  be  generalized  to 

Np  =  Z  Npk  (6.3) 

where 

NPk  ’  "k *  W  t6-4> 
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The  symbols  have  the  same  meaning  as  the  corresponding  symbols  without 
the  subscript  k,  except  that  they  now  represent  the  variables  for  a 
specific  nuclei  size  category  (which  is  characterized  by  a  mean 
radius). 

It  follows  that  if  Equation  (6,3)  is  an  adequate  scaling  parameter 
for  cavitation  inception  in  jets,  then  cr^  would  be  a  function  of  Np 
only.  That  is  to  say  a  universal  curve  would  be  obtained  when  a. 
is  plotted  against  Np,  irrespective  of  flow  conditions.  The  results  of 
the  present  work  were  used  to  ascertain  if  Np  is  a  useful  scaling 
parameter. 

The  values  of  Np  corresponding  to  the  various  measured  a.  's 
were  computed  at  jet  velocities  of  10  m/s,  20  m/s  and  30  m/s.  The 
method  by  which  these  values  were  obtained  is  shown  in  Appendix  2 
and  only  the  results  are  discussed  here.  However,  it  should  be  pointed 
out  that  two  assumptions  were  made  to  enable  these  calculations  to 
be  carried  out  and  they  were: 

(i)  that  the  peak  pressure  fluctuation  intensities  scaled  with 
the  dynamic  head  of  the  jet  at  exit;  and 
(ii)  that  the  probability  distribution  curve  of  the  peak  pressure 
fluctuations  was  the  same  for  all  the  three  jets. 

These  assumptions  had  to  be  made  because  the  pressure  measurements 
were  performed  at  only  one  velocity  (=  13  m/s)  and  furthermore  no  data 
were  available  on  the  largest  jet  for  reasons  already  discussed  in  the 
previous  chapter. 


The  calculated  results  are  summarized  graphically  in  Figures  6.1 
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to  6.3.  It  can  be  inferred  from  these  plots  that  in  the  range  of 
velocities  of  10  m/s  to  30  m/s,  the  measured  a^‘s  do  scale  with  Np  , 
for  a  constant  size  jet.  However,  when  the  data  on  the  three  jets 
were  plotted  on  the  same  graph,  the  data  points,  especially  those  for 
the  6.35  nm  jet,  do  not  all  fall  on  the  same  curve.  It  is  conceivable 
that  the  probability  distribution  curve  of  the  peak  pressure  fluctua¬ 
tions  is  not  universal  and  is  a  function  of  the  jet  size  or  some  yet 
undefined  parameter.  Consequently,  this  would  mean  that  the  computed 
Np's  for  the  largest  jet  were  not  scaled  in  the  right  proportion.  It 
would  be  premature  to  say  that  this  is  indeed  the  reason  for  the  scat¬ 
ter  in  the  data  points  without  further  investigation.  Nevertheless, 
there  is  a  strong  indication  that  a..  is  a  function  of  Np. 

We  have  already  mentioned  in  Section  2.6  that  in  the  current 
investigation,  inception  was  called  when  5  cavitation  events  or  more 
per  second  were  detected.  It  was  suggested  that  it  will  be  more  ap¬ 
propriate  to  scale  the  count  of  events  to  the  volume  in  the  cavita- 
ting  region  of  the  jet.  Thus  in  a  larger  jet,  inception  will  be  called 
at  a  correspondingly  higher  count.  This  means  that  in  the  present 
case,  cavitation  inception  in  the  largest  jet  should  only  be  called 
when  about  50  cavitation  events  or  more  are  measured.  The  question 
at  this  point  is  whether  the  measured  c^'s  are  significantly  changed 
as  a  result  of  this  new  definition  of  the  Inception  point.  We  claim 
that  the  values  of  's  are  not  greatly  affected  because  it  has  been 
observed  that  once  cavitation  inception  occurs,  the  number  of  cavitation 
events  increases  exponentially  thus  implying  that  the  mean  static  pres- 
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sures  between  the  two  inception  points  (and  consequently  were  not 
that  different.  This  is  especially  true  for  jets  with  a  high  nuclei 
population,  A  statistical  analysis  will  now  be  performed  to  verify 
the  above  claim. 

It  is  to  be  expected  that  the  probability  of  the  occurrence  of 
cavitation  events  is  governed  by  the  number  of  nuclei  in  the  cavitation- 
prone  region  and  the  cumulative  probability  of  favorable  fluctuations, 
that  is  the  parameter,  Np,  as  discussed  earlier.  Since  the  occurrence 
of  these  events  is  a  Poisson  process,  the  probability  of  say,  m  events 
occurring  per  second  can  be  Written  as 


and 


PM 


m 

P(>m)  *  1  -  E  P(4>) 


(6.5) 


(6.6) 


where  P(*m)  denotes  the  probability  of  the  occurrence  of  m  events  or 
more  per  second.  As  a  conservative  estimate  we  will  take  Np  *  500. 
Therefore,  from  Equations  (6.5)  and  (6.6)  we  have 

P{>5)  -  1  -  e'500  [500  ♦  |“£  ♦  §?°L  ♦  |°SL] 
or 


} 


P(>5)  *  1 
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because  the  product  of  the  term  in  the  bracket  and  the  exponential  term 
is  very  small  in  comparison  to  unity.  Similarly, 


PC>501  *  1  -  e 


*500 
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5Q0  + 


50Q2  , 

TT“  + 
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=  1 

As  demonstrated  above,  the  occurrence  per  second  of  these  cavitation 
events  has  very  nearly  the  same  probabilities  thereby  supporting  the 
claim  that  the  measured  inception  indices  in  the  present  investigation 
were  not  significantly  affected  by  the  technique  chosen.  However, 
it  can  also  be  deduced  from  the  above  analysis  that  if  Np  is  small  in 
comparison  to  the  number  arbitrary  assigned  as  the  inception  point,  then, 
the  difference  in  value  of  the  measured  by  virtue  of  this  different 
definition  of  the  Inception  point  would  be  significant. 

Another  interesting  finding  of  the  current  studies  is  that  the 
cavitation  inception  number  for  a  particular  nozzle  size  is  independent 
of  the  jet  exit  Reynolds  number,  at  least  over  the  range  of  velocities 
investigated.  A  possible  reason  for  this  observation  is  discussed  in 
the  next  paragraph. 

Consider  the  case  of  a  jet  of  an  arbitrary  size.  As  mentioned 
before,  inception  would  only  occur  when 


or 
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pn  i  PV  '  3* 

Tf  Inception  occurred  at  a  mean  pressure  denoted  by  PT,  then  the  minimum 
negative  pressure  fluctuation  required  to  initiate  cavitation  is  given 
by 


P*  -  CP*  -  PT1 


(6.7) 


Substituting  for  P£  in  Equation  (6.7)  gives 


p'  '  tPV  -  3R-  -  PT>  t6-81 

Furthermore  if  P'  *  KHp  V^2,  where  1C  is  a  constant,  then  Equation  C6.8) 
can  be  rewritten  as 


K  * 


(Py  PT1 

- 2 — 

*5P  V/ 


8S 


3p  RVh 


or 


ai 


(Ka  -  A0f) 


(6.9) 


where 
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PT-PV 

HP  V.2 
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3p  RV. 


and 


K  *  -K  *  constant 
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For  the  purpose  of  discussion,  Equation  6,9.  was  plotted  over  the  range 
of  nuclei  sizes  found  In  the  present  jets,  This  Is  shown  in  Figure 
6.5  for  an  arbitrary  value  of  K  ,  In  the  present  work,  ct.  was  measured 

a  1 

In  the  velocity  range  of  approximately  20  m/s  to  35  m/s.  It  can  be 
seen  from  this  figure  that  over  this  range  of  velocities,  ck  is  only  a 
weak  function  of  the  jet  exit  velocity,  and  the  maximum  change  in  the 
value  of  Act.,  is  0,05,  This  change  in  value  of  Act.  with  velocity  is  of 
the  same  magnitude  as  the  scatter  in  the  measured  ct. 's.  It  is  probably 
because  of  this  reason  that  the  measured  cr. 's  give  the  appearance  of 
being  independent  of  Reynolds  number  at  a  constant  jet  size. 

Since  K  is  arbitrary,  it  is  applicable  to  any  jet.  If  K,  had 

a  a 

been  independent  of  jet  size,  then  Figure  6.5  says  that  the  measured  ct^'s 
for  the  three  jets  investigated  would  have  the  same  ct..  's  since  they 
were  operated  over  the  same  velocity  range  and  had  similar  nuclei  size 
distribution.  However,  the  measured  ct^'s  are  otherwise  thus  suggesting 
that  K  could  possibly  be  a  function  of  jet  size  and  maybe  even  some 

a 

other  parameters.  It  should  be  emphasized  again  that  all  the  above 
arguments  are  only  valid  if  K  is  a  function  of  the  dynamic  head.  This 
information  is  presently  unavailable,  bearing  in  mind  that  K  is  a 

Q 

"scaling  constant"  for  the  peak  and  not  the  rms  pressure  fluctuations. 
Further  investigations  are  required  before  any  definite  conclusions  could 
be  drawn. 

Another  puzzling  observation  was  the  location  of  the  cavities  at 
inception.  As  has  been  pointed  out  in  the  last  chapter,  in  the  two 
smaller  jets,  the  cavities  were  first  seen  in  the  region  beyond  the 


potential  core;  whereas,  in  the  largest  jet  the  cavities  were  just  as 
likely  to  occur  in  the  shear  layer  wedged  Between  the  potential  core 
and  the  entrainment  zone.  As  a  result  of  this  observation  ,  we 
believe  that  it  is  appropriate  for  the  cavitation-prone  region  of  any 
jet  to  be  divided  into  two  sections;  namely,  the  shear  layer  region  and 
the  region  beyond  the  potential  core  which  is  referred  to  as  the  transi¬ 
tion  region  of  the  jet  in  some  literature.  Furthermore,  we  would 
theorize  that  the  factor  that  determines  where  the  cavity  is  first  seen 
is  the  ratio  of  the  Np  in  the  shear  layer  to  that  in  the  transition 
zone.  It  is  conceivable  that  when  the  ratio  is  less  than  a  certain  yet 
undetermined  value,  inception  would  always  occur  in  the  transition  zone 
of  the  jet.  On  the  other  hand,  this  ratio  of  the  probable 
cavitation  occurrence  parameter,  Np  ,  could  be  so  large  that  the  cavities 
at  inception  always  appear  in  the  shear  layer,  for  all  flow  conditions. 
Between  these  two  extremes,  there  is  probably  an  intermediate  range 
where  cavities  could  appear  in  either  or  both  regions  in  the  jet. 
Unfortunately  a  quantitative  analysis  similar  to  that  performed  earlier 
could  not  be  carried  out  here  due  to  a  lack  of  information  on  the 
pressure  fluctuation  intensities  in  the  shear  layers  of  the  present 
jets.  However,  if  we  were  to  use  the  size  of  the  jet  as  a  rough 
guide-line,  the  present  results  would  suggest  that  for  jets  smaller 
than  4.76  mm,  the  cavities  would  almost  always  occur  In  the  transition 
zone.  In  jets  larger  than  4.76  mm,  the  cavities  are  just  as  likely 
to  appear  In  the  shear  layer.  The  upper  bound  of  this  intermediate 
range  is  probably  defined  by  a  jet  of  25.4  mm  size  where  it  has  been 
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reported  by  Arndt  u  L  for  example  that  the  cavities  were  seen  in  the 
shear  layer.  Again  there  is  a  need  for  further  work  in  this  area 
before  the  actual  cause  of  this  phenomenon  is  established. 

The  results  of  the  present  studies  as  well  as  those  of  others  have 
shown  that  at  a  high  dissolved  air  content,  inception  in  the  jet  oc¬ 
curred  at  a  relatively  high  mean  pressure  in  the  jet.  In  such  a 
situation,  only  a  few  locations  in  the  jet  experienced  local  pressure 
fluctuations  that  were  intense  enough  to  cause  cavitation.  Consequently, 
the  cavities  at  inception  only  occurred  in  these  few  positions.  How¬ 
ever,  at  a  low  air  content,  the  mean  pressure  at  inception  was  con¬ 
siderably  lower  which  resulted  in  many  locations  having  local  pressures 
that  were  conducive  to  the  explosive  growth  of  the  nuclei.  Thus,  the  I 

cavities  were  much  more  extensive  than  the  former  case.  These  ex¬ 
planations  are  consistent  with  the  results  of  our  flow  visualization. 

Up  to  now,  we  have  only  presented  a  discussion  on  cavitation  in-  , 

ception  in  the  jets  where  the  cavitation  nuclei  are  not  generated  by 
any  external  sources  but  are  from  the  dissolved  air  content  in  the  water. 

A  discussion  on  the  behavior  on  the  same  jets  when  flow  seeding  was 

performed  will  now  be  considered,  j 

6.3  Seeded  Jets 

It  is  currently  unknown  if  the  pressure  field  In  the  jet  was 
affected  by  the  introduction  of  electrolysis  bubbles,  especially  in 
cases  when  the  bubbles  were  present  in  large  numbers.  However,  for  the 
sake  of  this  discussion,  it  will  be  assumed  that  the  pressure  field 
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remains  the  same  as  in  the  unseeded  jet,  Thus  it  is  presumed  that 
the  only*  effect  of  the  seeding  was  to  alter  the  nuclei  population. 

We  have  already  pointed  out  that  when  flow  seeding  was  performed, 
the  nuclei  were  so  abundant  that  the  counting  of  these  nuclei  from  the 
holograms  was  impractical.  Nevertheless  it  is  still  desirable  to  have 
some  estimate  on  these  nuclei  populations.  From  Kuiper's  work 
the  volume  of  hydrogen  gas  produced  per  unit  length  by  an  electrolysis 
current  at  a  water  temperature  of  23°  C  is  given  by 

G  =  1 ,27xl0"4  jip-  m3/s  £6.10) 

(d  g 

where  I  *  current  through  the  wire  in  amperes 
=  the  length  of  the  wire  in  meters,  and 
Pg  =  the  gas  pressure  in  N/m^. 

The  gas  pressure  Pg  can  be  found  from  the  static  equilibrium  equation, 
viz 

Pg  '  (p  -  PV  -  ir> 

where  P  is  the  static  pressure  at  the  position  of  the  wire.  It  is 
surprising  to  note  from  Equation  (6.10)  that  the  volume  of  gas  produc¬ 
tion  is  independent  of  the  diameter  of  the  electrode.  To  obtain  an 
order  of  magnitude  estimate  on  the  nuclei  population,  we  shall  further 
assume  that  all  the  nuclei  are  of  the  same  size  and  are  10  ym  in 
diameter.  Hence,  If  n^  denotes  the  number  of  bubbles  produced  per 
unit  time  by  the  current  I,  then, 
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Substituting  for  6  in  Equation  (6,111  and  rearranging  we  have 

n  »  2.43xlQU  (6.12a) 

y  \ 

But  as  mentioned  in  previous  chapters,  =  0,5  mm.  As  can  be  infer¬ 
red  from  Equation  (6. 12a), the  number  of  bubbles  generated  is  dependent 
on  the  gas  pressure  which  in  turn  is  determined  by  the  velocity  of  the 
jet.  However,  for  convenience,  n^  was  computed  at  only  one  represen¬ 
tative  velocity  of  40  ft/s  (13  m/s).  This  corresponded  to  a  gas 
pressure  of  approximately  20  psia  0*35  atm,).  Thus  Equation  (6.12a) 
becomes 


n.g  *  3. 52x1 09  I 


(6.12b) 


Hence  it  follows  that  the  total  number  of  nuclei  that  are  present  in 
the  cavitation-prone  region  at  any  instant  is  given  by 


n 


nd 


+ 


where  n^  =  number  of  nuclei  normally  present  in  the  jet  without  flow 
seeding 


Q  *  volume  of  the  cavitation-prone  region 
and  Qj  *  volume  efflux  of  the  nozzle  per  second. 

Proceeding  In  a  similar  manner  as  has  already  been  discussed  for 
the  unseeded  jets,  the  measured  incipient  cavitation  indices  at  various 
magnitudes  of  the  electrolysis  current  were  plotted  against  the  para- 
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meter  Np  for  all  the  three  jets  concerned,  This  Is  shown  graphically 
In  Figure  (6.6),  Again,  It  can  Be  seen  that  while  the  data  for  the 
two  smaller  jets  lie  on  the  same  carve,  those  for  the  largest  jet  do 
not.  The  probable  explanation  for  this  discrepancy  has  already  been 
discussed  in  the  previous  section. 

It  can  be  easily  shown  that  when  flow  seeding  was  performed,  the 
number  of  bubbles  introduced  into  the  flow  were  far  more  numerous  than 
those  that  are  normally  present  in  the  jet.  Consequently  one  would 
not  expect  the  dissolved  air  content  of  the  water  to  have  any  signifi¬ 
cant  effect  on  the  cavitation  inception  number  under  these  circum¬ 
stances.  This  was  the  general  observed  behavior  of  the  jets  under 
investigation  as  was  seen  from  Figures  5,20  to  5.43.  There  were  some 
exceptions  where  the  dissolved  air  content  had  an  observable  effect 
on  the  measured  a. 's  although  no  regular  pattern  could  be  inferred. 

The  reason  for  this  discrepancy  in  the  data  is  not  known  at  present. 
Another  observation  of  the  present  results  is  that  do  not  scale 
with  the  jet  exit  Reynolds  although  in  some  instances  the  data  of  the 
jets  do  supplement  one  another. 

With  flow  seeding,  the  jet  cavitated  a  lot  sooner  and  as  a  result 
the  tests  had  to  be  confined  over  a  lower  range  of  velocity  which 
was  10  m/s  to  20  m/s.  Looking  back  at  Figure  (6.5)  it  can  be  seen 
that  over  this  range  of  velocity,  there  is  a  strong  dependence  of  ai 
on  the  jet  velocity  (or  Reynolds  number  at  a  constant  jet  size). 
Therefore  we  would  expect  this  dependence  to  show  up  despite  the 
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scatter  of  the  data  points,  which,  is  the  case  as  can  he  seen  in  Figures 
5.8(al  to  5,10(41,  It  should  be  stressed  at  this  point  that  the 
above  explanation  is  valid  only  if  the  assumptions  made  for  the  plot 
of  Figure  6,6  are  true.  The  experimental  results  also  showed  that  for 
a  constant  jet  size  and  at  a  fixed  velocity,  a.,  increases  with  increas¬ 
ing  electrolysis  current.  The  reason  for  this  behavior  can  be  ex¬ 
plained  as  follows.  When  the  flow  was  seeded  with  more  and  more 
nuclei,  there  was  an  increasing  tendency  for  cavitation  inception 
to  be  caused  by  higher  and  higher  fluctuation  intensities  although 
the  frequency  of  occurrence  of  these  intensities  were  less  than  the 
lower  ones.  Consequently  increasing  values  of  ai  were  recorded. 

In  addition,  a  "saturation"  point  was  reached  when  the  nuclei  were  so 
numerous  that  inception  in  the  jet  occurred  at  the  very  first  oppor¬ 
tunity.  Beyond  this  point,  the  measured  a,  did  not  exhibit  any  sig¬ 
nificant  change  in  value  with  further  increase  in  the  number  of  electro¬ 
lysis  bubbles. 

We  have,  up  to  this  moment,  only  discussed  the  results  of  the 
present  investigations  with  little  to  no  references  made  to  the  work 
of  others.  Our  attention  will  now  be  turned  to  this  area. 

6.4  Comparison  of  Experimental  Results 
6.4.1  Nuclei  Measurements 

It  has  already  been  stated  earlier  that  so  far  measurements  of 
the  distribution  of  the  nuclei  population  in  submerged  jets  have  not 
been  available.  However,  numberous  investigators,  for  example.  Gates 
and  Acosta (46),  Peterson  et  al  and  Keller  and  Weitendorf  (6®)  have 
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made  nuclei  measurements  in  water  tunnels  under  various  flow  conditions. 
It  can  be  seen  by  comparing  Figure  6.7  and  Figure  6.8  that  the  results 
of  the  current  work  compare  favorably  with  those  of  the  others  except 
that  large  size  nuclei  were  not  found  in  our  jets.  This  is  because 
the  present  measurements  were  made  in  a  jet  discharging  into  an  other¬ 
wise  stationary  medium;  whereas,  the  other  measurements  were  performed 
in  flowing  water  tunnels  which  were  at  times  under  cavitating  condi¬ 
tions. 

In  the  preliminary  work  by  Ooi  and  Katz  it  was  reported 

that  the  size  and  distribution  of  the  nuclei  population  were  dependent 
on  the  cavitation  number  which  is  contrary  to  the  present  observa¬ 
tion.  Since  the  measurements  were  made  on  the  same  test  facility, 
we  believe  that  this  discrepancy  in  the  results  is  due  to  the  dif¬ 
ferent  test  procedures.  Recall  that  in  the  current  investigation, 
the  pressure  in  the  tunnel  was  always  raised  to  its  initial  value 
and  the  tunnel  was  permitted  to  sit  at  that  pressure  for  some  time 
before  the  next  test  run  was  made.  This  procedure  was  not  adopted  in 
the  earlier  reported  work.  It  is  thus  obvious  that  there  is  a  need 
for  all  investigators  to  monitor  the  nuclei  population  distribution 
in  their  jets  especially  if  different  test  procedures  are  followed. 

6.4.2  Photographic  Investigation  of  the  Non-Cavltating  Jets  - 
Coherent  Structures 

The  existence  of  large  quasi  two-dimensional  coherent  structures 
in  shear  flows  was  first  reported  by  Brown  and  Roshko  Since 
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then,  numerous  other  investigators,  for  example,  Winant  and  Browand^^, 
Konrad^73}  Davies  and  Yule^^and  Perry  and  Watmuff^^have  directly  or  in¬ 
directly  established  the  presence  of  these  structures  in  shear  layers, 
wakes  and  jets.  However,  Chandrsuda  et  al^  ^7^  and  Castro  and 
Bradshaw  suggested  that  these  large  scale  eddies  were  neither 
quasi-deterministic  nor  two  dimensional.  Both  groups  of  investigators 
claimed  that  the  observed  features  of  Brown  and  Roshko  were  "the 
notoriously  prolonged  effects  of  transition"  of  initially  laminar 
layers.  Chandrsuda  et  further  suggested  that  the  flow  would 
become  highly  three-dimensional  whenever  the  free  stream  turbulence 
was  high.  However,  these  claims  have  been  disproved  by  Wygnanski 
et  al^  who  concluded  that  "the  two-dimensional  character  of  the 

coherent  eddies  perseveres  in  spite  of  strong  external  small-scale 
buffeting,  and  the  Brown-Roshko  structure  may  be  more  common  in 
practice  than  was  hitherto  observed". 

In  the  present  study,  these  large  structures  could  only  be  ob¬ 
served  up  to  a  relatively  low  Reynolds  number  (based  on  axial  posi¬ 
tion)  of  approximately  45,000.  An  example  of  these  structures  is 
shown  in  Figure  6.9.  Unlike  the  structures  in  the  plane  shear  layer 
described  by  Brown  and  Roshko,  the  present  large  scale  structures  only 
remained  distinct  up  to  the  second  or  third  vortex.  Beyond  that 
the  flow  became  three-dimensional  with  a  multitude  of  small  grainy 
structures. 


t 


As  was  pointed  out  in  the  last  chapter,  cellular  structures  were 
seen  just  downstream  of  the  laminar  region  of  the  jet.  These  struc- 
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tures  were  also  seen  by  Konrad  v  .  It  is  believed  that  these 
cellular  structures  begin  as  hairpin  vortices  and  as  they  proceed 
further  downstream,  they  become  unstable  and  produce  the  three- 
dimensional  small  scale  structures  in  the  jet. 

The  work  of  Brown  and  Roshko  has  established  that  the  longitu¬ 
dinal  increase  in  the  scales  of  the  large  eddy  structures  and  con¬ 
sequently  the  spreading  rate  of  the  shear  layer  was  accomplished  by 
the  coalescence  of  two  or  more  of  the  neighboring  coherent  structures. 
Based  on  their  schlieren  photographs.  Brown  and  Roshko  estimated  that 
the  slope  of  the  spread  of  the  shear  layer  was  0.38.  In  the  present 
work,  slopes  of  0.33,  0.24  and  0.43  were  obtained  for  the  3.17  mm, 

4.76  mm,  and  6.35  mm  jets,  respectively.  The  reason  for  the  different 
spread  rates  is  presently  unclear.  It  has  been  suggested  that  the 
initial  condition  of  the  jet  could  be  a  determining  factor. 

6.4.3  Pressure  Fluctuation  Measurements 

All  previous  pressure  fluctuation  measurements  were  made  in  an 
Eulerian  frame  with  intrusive  pressure  transducers.  Although  such 
measurements  do  give  information  concerning  the  pressure  field  in 
the  jet,  they  are  not  directly  applicable  to  the  problem  of  cavita¬ 
tion  which  is  a  Lagrangian  problem  (i.e.,  the  fluctuations  moving 
with  the  fluid).  It  is  our  belief  that  the  bubble  injection  technique 
used  in  the  present  work  and  described  in  detail  in  Chapter  4  repre¬ 
sent  the  first  attempt  ever  in  measuring  the  pressure  fluctuations  in 
a  Lagrangian  frame.  The  present  method  has  the  added  advantage 
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that  the  "probes"  do  not  cause  a  significant  flow  disturbance  since 
they  are  a  part  of  the  flow.  Also  this  method  provides  a  means  where 
the  instantaneous  pressures  at  numerous  positions  in  the  jet  could  be 
simultaneously  measured,  which  is  not  possible  with  the  conventional 
method  by  virtue  of  the  fact  that  the  flow  would  be  drastically  changed 
by  the  inserted  pressure  transducers. 

It  has  already  been  pointed  out  in  the  course  of  this  chapter  that 
cavitation  is  initiated  when  the  local  pressure  dips  to  or  below  the 
vapor  pressure.  Thus  a  knowledge  of  the  peak  temporal  pressure  fluc¬ 
tuations  would  be  very  useful.  However,  with  the  exception  of  the 
present  investigation,  such  information  has  been  sorely  lacking.  All 
the  measurements  have  dealt  with  the  rms  of  the  pressure  fluctuations 
although  Arndt  and  Ippen  ^  and  Rouse  for  example  have  noted 

that  peak  fluctuations  could  be  as  high  as  7  times  the  rms  level  in  a 
turbulent  boundary  layer  and  10  times  ip  a  jet.  In  the  present  in¬ 
vestigation,  the  maximum  peak  fluctuation  intensities  are  only 
about  2  to  3  times  that  of  the  computed  rms  values.  Various  possible 
reasons  for  this  discrepancy  come  to  mind.  The  fact  that  the  present 
measurements  were  made  in  water  whereas  all  the  other  measurements 
were  in  air  may  account  for  some  of  this  difference.  It  has  been  dis¬ 
cussed  in  earlier  sections  that  the  intensity  of  the  pressure  fluc¬ 
tuations  could  very  well  be  a  function  of  the  jet  size  and  velocity. 
Both  these  factors  could  be  additional  reasons  for  the  discrepancy 
in  the  results.  In  any  event,  the  peak  pressure  fluctuation  inten¬ 
sities  in  the  present  jets  are  high  enough  to  initiate  cavitation 
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inception  if  a  nucleus  is  present  at  the  right  instant. 

To  have  an  idea  on  the  rms  fluctuation,  the  jet  was  divided  into 
various  regions  as  shown  in  Table  6.1.  The  calculation  was  then  per¬ 
formed  using  the  expression: 
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where  < - *■>  =  rms  of  the  dimensionless  pressure  fluctuation  in  a 

*P  vj 

particular  region  of  the  jet. 

m  *  total  number  of  data  points  in  the  same  region  of  the 

jet. 


P*  *  instantaneous  pressure  fluctuation. 

The  computed  results  are  tabulated  in  Tables  6.1  and  6.2  for  the 
polytropic  and  isothermal  cases,  respectively.  It  can  be  seen  from 
these  tables  that  there  are  two  entries  in  each  region  of  the  jet. 

The  first  entry  is  the  calculated  rms  fluctuation  when  all  the  measured 
data  were  used.  The  second  entry  corresponds  to  the  calculated  rms 
value  when  a  few  of  the  data  points  were  not  included  in  the  calcula¬ 
tion  because  they  exhibited  values  that  were  significantly  different 
from  the  majority  of  the  data  in  that  group.  A  quick  comparison  of 
these  two  entries  clearly  show  that  in  many  instances,  there  is  a 
sizeable  difference  in  the  two  rms  values.  This  implies  that  the 
calculated  values  are  not  too  accurate  because  of  the  limited  number 
of  data  points.  It  should  also  be  pointed  out  that  the  numbers  in 


0<r/D<0.5  0.5<r/D<1.0  1.0<r/D<1.5_  l,5<r/D<2,0  2,Q<r/D<2,5  2,5^r/D<3,0 


Table  6.1  A  tabulation  of  the  rms  pressure  fluctuations  In  various  regions  of  the  present  jet. 
(poly tropic  Index  =4/3) 


Qir/D<P.5  Q,5£r/P<1.0  1.0<r/D<1.5  l.S<r/D<2,0  2.0<r/D<2.5  2.5<r/D<3.0 

0.19(95)  0.33(38)  0.49(13) 


Table  6.2  A  tabulation  of  the  mns  pressure  fluctuations  In  various  regions  of  the  present  jet. 
(Isothermal  case) 


parentheses  in  the  two  tables  represent  the  sampling  size. 

In  the  rms  pressure  fluctuation  measurements  made  to  date, 
the  results  show  that  the  magnitude  of  luctuations  decreases  with 
increasing  distance  from  the  jet  exit.  However,  this  trend  is  not 
apparent  in  the  present  work.  The  reason  for  this  discrepancy  is  not 
known  at  present.  In  addition  the  current  rms  level  is  considerably 
higher  than  those  reported  by  Sami  et  al_  and  Barefoot 
Their  lower  rms  measurements  could  be  due  to  the  effect  of  spatial 
averaging  since  the  transducers  they  used  were  large  by  comparison. 
Needless  to  say,  further  work  has  to  be  undertaken  before  some  of 
these  differences  in  experimental  findings  can  be  explained. 

6.4.4  Cavitation  Results 

A  careful  survey  on  the  published  work  in  jet  cavitation  revealed 
that  this  area  is  currently  in  a  very  disorganized  state.  Despite  the 
fact  that  measurements  were  usually  performed  using  different  experi¬ 
mental  procedures  under  various  flow  configurations,  many  investigators 
have  not  carefully  taken  into  account  the  flow  environment  and  the 
flow  dynamics  associated  with  the  jets  they  were  working  on.  For 
example,  only  a  handful  of  workers  have  carried  out  flow  visualization 
studies  and  pressure  fluctuation  measurements  in  jets.  Nuclei  popula¬ 
tion  density  measurements  have  also  been  seriously  lacking.  As  a 
result  of  this  oversight,  it  is  often  difficult  if  not  impossible  to 
compare  the  results  of  the  various  investigators  especially  since 
o^'s  were  usually  measured  in  different  sized  nozzles  and  at  various 
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dissolved  air  content.  Nevertheless  we  shall  present  a  very  brief 
discussion  on  the  comparison  of  the  present  results  with  the  work  of 
others  below.  This  discussion  will  be  limited  to  jet  cavitation  only 
since  other  shear  flow  results  can  be  found  in  Arndt  and  George 
and  Arndt 

The  incipient  cavitation  number  for  a  3.17  mm  and  a  6.35  mm  jet 
was  reported  to  be  0.25  and  0.55,  respectively,  by  Lienhard  and 
Stephenson  ^8).  Although  they  did  not  measure  the  air  content  in 
their  jet,  we  believe  the  jet  could  be  at  most  near  saturation.  This 
means  that  their  a_.'s  are  considerably  higher  than  the  presently 
measured  values  of  0.16  and  0.27  for  the  unseeded  jet  at  14  ppm. 

The  higher  values  of  a.-’s  could  be  due  to  the  following  reasons. 

It  has  been  shown  from  schlieren  photographs  that  the  present 
jets  are  laminar  at  exit.  On  the  other  hand  it  is  likely  that  the 
jets  of  Lienhard  and  Stephenson  were  turbulent.  If  this  is  true 
then  the  higher  could  be  attributed  to  the  higher  nuclei  population 
in  the  cavitation-prone  region  of  the  jet  since  Hussain  and  Zedan 
have  demonstrated  that  the  spread  rate  of  an  initially  turbulent  jet 
is  greater  than  a  laminar  one.  The  difference  in  the  results  could 
also  be  caused  by  the  different  intensity  of  pressure  fluctuations 
since  it  is  conceivable  that  they  are  governed  by  the  initial  con¬ 
ditions  of  the  jet. 

Another  disagreement  in  the  results  is  the  location  of  the 


cavities  of  inception  in  the  jet.  It  has  either  been  stated  or  im¬ 
plied  in  a  number  of  sources  (e.g.,  Arndt  and  Arndt  and  George^ 10 b 
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that  the  cavities  at  inception  are  first  seen  in  shear  layer  of 
the  jet,  which  is  not  always  the  case  in  the  present  jets.  The 
probable  explanations  for  this  discrepancy  have  already  been  discussed 
earlier  and  will  not  be  repeated  here.  However,  we  would  like  to 
comment  on  Sami  et  al_'s  assumption  that  the  cavities  usually 

occur  at  the  cores  of  the  turbulent  eddies.  Our  extensive  flow 
visualization  studies  in  the  present  jets  have  not  given  such  an  indi¬ 
cation.  Of  the  hundreds  of  schlieren  photographs  that  were  taken 
at  inception,  only  one  showed  the  occurrence  of  a  cavity  in  the  core 
of  a  vortex.  This  is  shown  in  Figure  6.10.  Finally,  it  should  be 
also  pointed  out  that  contrary  to  the  findings  of  Baker  et  al^  (I4), 
the  measured  cavitation  index  of  the  present  jet  was  a  function  of 
the  dissolved  air  content  even  when  the  jet  was  undersaturated.  The 
reason  for  this  discrepancy  is  still  not  known. 

6.5  Proposals  for  Future  Work 

One  of  the  assumptions  made  in  the  course  of  the  current  dis¬ 
cussion  is  that  the  probability  distribution  curve  of  the  peak  pres¬ 
sure  fluctuation  intensities  is  universal.  Questions  were  raised 
as  to  the  validity  of  this  assumption  since  the  distribution  curve 
could  very  well  be  a  function  of  the  jet  size,  the  jet  velocity  and 
the  exit  condition.  All  of  these  questions  could  be  answered  by 
performing  a  series  of  pressure  measurements  on  jets  of  different 
sizes  over  a  range  of  velocities  and  turbulence  level  at  exit.  Recall 
also  that  the  present  pressure  measurements  were  confined 


to  the  transition  region  of  the  jet.  Obviously  these  measurements 
could  not  be  accurately  used  to  predict  cavitation  occurrences  in  the 
shear  layer  since  the  distribution  curve  could  be  different.  Conse¬ 
quently,  it  would  be  advantageous  to  extend  these  measurements  to  this 
region  of  the  jet. 

It  was  also  postulated  that  the  factor  that  determines  where  the 
cavity  is  first  seen  at  cavitation  inception  is  the  ratio  of  the  Np 
in  the  shear  layer  to  that  in  the  transition  zone  of  the  jet.  This 
could  be  established  one  way  or  another  by  using  the  results  of  the 
above  proposed  studies. 


PROBABLE  CAVITATION  OCCURRENCE- 


Fig.  6.1.  A  plot  of  the  probable  cavitation  occurrence  parameter,  Np, 

versus  <j,  at  different  jet  exit  velocities  for  the  3.17mm  jet. 


PROBABLE  CAVITATION  OCCURRENCE 
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Fia  62  A  plot  of  the  probable  cavitation  occurrence  parameter,  Np, 

versus  ®1  at  different  jet  exit  velocities  for  the  4.76mm  jet. 
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PROBABLE  CAVITATION  OCCURRENCE- 
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for  the  equation  ■  Ka  -  (8S/3pRVj2). 


PROBABLE  CAVITATION  OCCURRENCE- 
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CAVITATION  INCEPTION  NUMBER-^ 


Fig.  6.6.  A  plot  of  the  probable  cavitation  occurrence  parameter,  N_, 
versus  o,  at  a  jet  exit  velocity  of  20m/s  for  the  three  K 
jets  when 'flow  seeding  was  carried  out. 
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Fig.  6.7.  A  plot  showing  the  nuclei  distribution  from  various  sources 
[Katz  and  Acosta  (1981)]. 
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MEAN  NUCLEI  RADIUS  -  R  (^.m  ) 

Fig.  6.8.  A  plot  showing  the  nuclei  distribution  In  the  present  jets  at 
different  dissolved  air  contents. 


SUMMARY  AND  CONCLUSIONS 


A  new  technique  employing  specially  tailored  air  bubbles  has 
been  successfully  used  to  map  out  the  temporal  pressure  field  in  the 
jet.  The  results  showed  that  negative  peak  fluctuation  of  as  low  as 
120  percent,  and  positive  peak  of  as  high  as  160  percent  of  the  dynamic 
head  existed  in  the  jet.  The  fluctuation  intensities  were  found  to 
have  a  slightly  skewed  bell  shape  distribution. 

Flow  visualization  revealed  that  cavitation  inception  occurred 
with  a  sudden  burst  of  cavities  within  the  jet.  The  location  of 
these  cavities  were  dependent  on  the  jet  size  but  were  never  observed 
beyond  16  jet  diameters  downstream  of  the  nozzle  lip.  In  general, 
these  cavities  did  not  start  in  the  cores  of  the  vortices. 

The  cavitation  inception  numbers  in  all  the  jets  were  found  to  be 
independent  of  the  velocity  when  the  jets  were  not  seeded.  The  dis¬ 
solved  air  content  had  a  definite  influence  on  the  inception  number. 

With  flow  seeding,  the  inception  number  exhibited,  in  most  cases,a 
linear  dependence  on  the  exit  velocity.  This  dependency  on  velocity 
was  found  to  increase  with  increased  nuclei  population.  The  present 
results  also  indicated  that  for  a  fixed  jet  size,  is  a  function 
of  a  newly  defined  paramet  Np,  called  the  probable  cavitation 
occurrence  parameter.  However,  we  were  unable  to  determine  if  this 
scaling  parameter  is  universal  from  the  present  experiments. 


A  survey  of  published  literature  showed  that  few  investigators 
have  actually  tried  to  correlate  their  cavitation  results  to  the  pres- 
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sure  field  in  the  jet.  In  addition,  many  of  these  workers  have 
neglected  to  make  careful  flow  visualization  studies  or  to  document 
the  nuclei  population  in  their  jets.  It  is  our  belief  that  until  and 
unless  such  studies  are  undertaken,  the  many  questions  that  exist 
will  never  be  satisfactorily  answered. 
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APPENDIX  1 

ORDER  OF  MAGNITUDE  ANALYSIS  OF  THE  TERMS  IN  THE  RAYLEIGH-PLESSET  EQUATION 
The  Rayleigh-Plesset  equation  is  given  by 


R  R 


♦  f  R2 


1  Cp  - 

p  L  c 


p»(t) 
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4vRt 

in 


(A.l ) 


If  the  circular  frequency  of  the  pressure  fluctuation  is  given  by  <*>, 
then  we  would  expect  the  change  In  R  to  be  proportional  to  u>  too.  Let 
R  be  given  by 


R  *  RQ  +  ee"Jwt  (A. 2) 
where  RQ  Is  the  initial  bubble  radius  and  e  is  the  amplitude  of  oscilla¬ 
tion  from  the  mean  radius.  It  should  be  pointed  put  that  e  could  be  of 
the  same  order  of  magnitude  as  RQ.  Therefore,  it  follows  from  Equation 
(A. 2)  that  q  r\j  ^ 


Furthermore,  we  define 


e  *  »  0(1 ) 

Ko 

w  *  — -  a  0(5),  6  «  l.(u)  is  the  natural  frequency  of  the 
o  0 

bubble).  Thus,  the  left-hand  terms  of  equation  (A.l)  become 
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In  addition 

pressure  term  ^  (u  ■  w  R  ) 

p  uo 
** 

viscous  term  ^  e  ui  *  0(5 )  , 

"R^ 

and  surface  tension  term  ^  i-  (VL  =  Weber  number  »  -^-7)  . 

We  e  pRu£ 

By  this  order  of  magnitude  arguments  we  can  conclude  that  if  the  time 
scale  of  the  pressure  fluctuations  is  small  in  comparison  with  the 
response  time  of  the  bubble,  then  the  Rayleigh-Plesset  equation  can  be 
simplified  to: 

0  •  pc  -  P.(t)  - 

or  _ 

p-<*>  '  Pc  -  “ 

*  Pg(T)  +  PV(T)  -  .  (A.3) 


176 


» 


» 


t 


♦ 


♦ 


APPENDIX  2.  COMPUTATION  OF  THE  PROBABLE  CAVITATION 
OCCURRENCE  PARAMETER,  Np 


A2.1  Objective 

The  aim  here  Is  to  determine  if  the  presently  measured  incipient 
cavitation  index  is  a  function  of  the  parameter  Np  as  defined  in 
Equation  (A2.1)  below.  If  it  turns  out  that  a  unique  curve  is  obtained 
when  c^.  is  plotted  against  Np,  then  this  would  strongly  suggest  that 
Equation  (A2.1)  is  an  accurate  description  of  the  scaling  law  that 
governs  cavitation  in  jets.  Before  the  plot  can  be  performed,  the 
values  of  Np  corresponding  to  the  various  measured  a- 's  for  the  three 
different  jets  operating  at  different  dissolved  air  contents, have 
to  be  computed.  This  was  performed  in  the  manner  shown  below  in 
Section  A2.2.  An  example  was  used  to  illustrate  the  steps  and  argu¬ 
ments  followed.  It  should  be  pointed  out  that  two  assumptions  were 
made  to  enable  the  calculations  to  be  carried  out.  They  were: 

i)  that  the  negative  peak  pressure  fluctuation  intensities 

2 

scaled  with  £p  IT  ;  and 

J 

i 1 )  the  probability  distribution  curve  of  pressure  fluctuation 
intensities  was  the  same  for  all  three  jets. 
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A2.2  Illustrative  Example 
Define 


Np  s  E  Npk  (A2.1 ) 

and 

NPk  “  nk  x  W  (A2.2) 

where 

nk  ■  total  number  of  nuclei  in  size  category  k  that  were  present 
in  the  cavitation-prone  region  of  the  jet;  and 

PjjPf)  a  cumulative  probability  of  favorable  negative  peak  pres¬ 
sure  fluctuation  intensities  that  would  cause  the  nuclei  in  the  said 
category  to  grow  explosively. 

Consider  the  case  of  the  3.17  mm  jet  at  a  dissolved  air  content 
of  14.1  ppm  and  at  a  velocity  of  20  m/s.  It  has  been  found  that 
inception  in  the  jet  occurred  at  a  mean  static  pressure, P-p of  5.04  psi 
(corresponds  to  ■  0.16).  In  addition,  nuclei  population  measure¬ 
ments  by  holography  revealed  that  nuclei  of  only  two  size  categories 
were  present  in  the  jets  investigated.  These  categories  corresponded 
to  mean  nuclei  radii  of  7.5  pm  and  1 5  van.  For  convenience,  subscripts 
1  and  2  will  be  used  In  Equation  (A 2.2)  to  denote  variables  that  are 
related  to  the  7.5  pm  and  15  pm  size  nuclei,  respectively. 

Bubble  dynamics  dictates  that  inception  would  only  occur  when 
the  Instantaneous  local  pressure  P£  satisfies  the  equation  below: 


(A2.3) 
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The  symbols  have  their  usual  meanings.  In  the  case  of  the  7.5  urn  size 
nuclei.  Equation  (A2.3)  is  satisfied  when 

P^s -1.44  psi 


while 


P£  i  "  0.52  psi  for  the  15  pm  size. 

The  negative  sign  indicates  that  the  fluid  is  under  tension.  Therefore, 

I 

the  minimum  pressure  f 1  uctuati on, P^, required  to  cause  the  onset  of 
cavitation  is  given  by 

Pmin  *  (pzmin  “  PT^ 

*  -6.48  psi  (for  the  7.5  pm  nuclei) 
and  Pmin  »  -5.56  psi  (for  the  15  pm  nuclei). 

In  this  example,  the  minimum  required  pressure  fluctuation  for  the  two 
size  categories  translates  into  22*  and  19%  of  the  dynamic  head, 
respectively.  From  the  probability  distribution  curve  in  Figure  5.45, 
the  cumulative  probabilities  are 

P^Pf)  *  0.33 
and 

P2(Pf)  *  0.43. 

For  convenience,  the  cumulative  probabilities  of  the  peak  fluctuation 
.  Intensities  from  negative  Infinity  to  specific  fluctuation  intensity 
levels  are  tabulated  as  shown  in  Table  A2.1  below. 


t 


p' 

P(Pf)xl02 

P' 

P(Pf)xl02 

ip  V2 

ip  V2 

o  > 

55.58 

-0.9  > 

4.93 

-0.1  > 

43.24 

-1.0  > 

3.08 

-0.2  > 

33.49 

-1.1  i 

1.85 

-0.3  > 

26.09 

-1.2  > 

0.86 

-0.4  > 

20.54 

-1.3  > 

0.43 

-0.5  > 

16.22 

-1.4  i 

0.18 

-0.6  > 

12.70 

-1.5  > 

0.06 

-0.7  > 

9.74 

-1.6  > 

0.00 

-0.8  > 

7.15 

le  A2.1.  A  tabulation  of  the  cumulative  probabilities  of  negative 
peak  pressure  fluctuation  intensities. 

In  addition, from  nuclei  measurements  we  have, 

n1  3  2.69xl03 
and 

n2  3  349  . 

Therefore,  Equations  (A2.1)  and  (A2.2)  imply  that 

Np  3  niXp-j (P^)  +  n2xP2(P^) 

3  (5690x0.33)  +  (349x0.43) 

3  1051.8 

Hence,  the  calculated  Np  corresponding  to  3  0.16  for  the  3.17  mm 
jet  is  1052. 


Proceeding  in  the  same  manner  illustrated  above,  the  values  of 
Np  for  the  other  measured  a^'s  were  generated.  The  calculations  were 
done  for  jet  velocities  of  10  m/s,  20  m/s  and  30  m/s.  The  calculated 
results  are  tabulated  as  shown  in  Tables  A2.2(a)  -  (c)  to  Tables 
A2.4(a)  -  (c)  below.  These  results  are  discussed  in  the  main  text. 


Cavitation  Air  Nuclei  size  =  7.5  pm  Nuclei  size  =  15  pm 

Inception  Control 


0.24  10.9  4.99  33.5  16.72  8.12  33.5  27.20  1944 
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